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ABSTRACT 

The  role  played  "by  the  f Tine tion- theoretic  Wiener-Hopf  techniques 
in  the  solution  of  certain  integral  equations  arising  from  the  application 
of  Green's  theorem  to  boundary  value  problems  le  well  known  and  has  recently 
been  extended.   In  the  present  work  we  discuss  and  illustrate  the  relevance 
and  applicability  of  these  techniques  when  the  method  of  separation  of  vari- 
ables in  cuiTvilinear  coordinates  is  employed.  For  the  same  boundary  various 
appropriate  or  semi-appropriate  coordinate  systems  may  be  employed,  while  on 
the  other  hand  the  Green's  function  integral  equation  is  invariant.   There 
is  seen  to  be  a  correspondence  between  the  two  methods,  which  are  expressions 
of  the  duality  between  Green's  function  and  eigenf unction  techniques.   Various 
aspects  of  this  correspondence,  as  well  as  the  effect  of  change  of  coordinete 
systems  ere  illustrated. 

The  resulting  relationships  shed  some  li^t  on  the  problem  of  recog- 
nizing the  relevancy  of  Wiener-Hopf  or  more  classical  techniques  to  certain 
integral  equations  and  serves  as  a  guide  in  the  discovery  of  new  problems  to 
which  the  method  (in  either  guise  mentioned  above)  may  be  applied.   Certain 
observations  with  regard  to  tiniqueness  also  result. 
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1,   Introduction 

The  exact  mathematical  theory  of  diffraction  has  had  a  history  dis- 
tinguished by  but  few  Bignlfieant  advances.  The  classical  method  of  separation 
of  variables  succeeds  readily  in  solving  the  problem  of  diffraction  of  light  or 
sound  waves  by  a  sphere.  Diffraction  by  other  obstacles  was  analyzed  in  an  ap- 
proximate manner,  but  as  late  as  the  close  of  the  nineteenth  century  Raylelgh  j_l] 
despaired  of  an  exact  solution  of  the  mathematical  problems  involved  in  dif- 
fraction by  a  str^iight  edge,  or  circular  aperture.     The  former  es- 
sentially two-dimensional  problem  was  first  solved  by  Sommerfeld  [2]  ,  by  an 
ingenious  extension  of  the  method  of  images  to  Rlemann  surfaces.  Apparently  this 
Bdthod  was  viewed  as  possessing,  to  a  too  great  extent,  the  character  of  an 
inimitable  tour  de  force,  for  the  problem  has  been  frequently  reworked  by  various 
means  of.  [3j,  [itj,  l5J  of  which,  however,  Bateman  [^5j  expresses  the  opinion 
that  they  all  require  special  insight.  One  feels  that  there  ought  to  exist  a 
unifying  method,  rather  than  a  disjoint  set  of  devices  each  applicable  but  once 
or  twice.   The  method  of  separation  of  variables  has  recently  been  applied 
to  diffraction  by  a  circular  disc  or  aperture  by  Bouwkamp[63  and  others  see 
for  example [21J .   It  had  already  been  often  applied  to  the  discussion  of  the 
diffraction  by  an  Infinite  ribbon,  or  by  a  slit  cf.  [^2   •   The  Mathieu  functions 
which  arise  have,  however,  been  in  an  Tinsatisfactory  state  until  lately.  During 
recent  years  various  sharp  improvements  in  technique  were  introduced,  among  them 
the  Hiener-Hopf  Method.   The  application  of  the  latter  to  diffraction  prob- 
lems appears  to  be  due  to  J.  Schwlnger  fsj,  in  the  United  States,  but  an  ex- 
ample of  it  also  occurs  In  work  by  E.  T.  Copson  {9J.  A  whole  class  of  electro- 
magnetic boundary  value  problems  Is  trested  thereby  in  a  novel  and  uniform 
manner.  A  generfilizatlon  due  to  G.  Carrier  fiol  will  be  discussed  below  in 

conjunction  with  this  matter. 

*  Numbers  in  square  brackets  refer  to  the  Bibliography. 
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This  procedure  of  Schwlnger  provided  the  initial  stimuluB  to  the 
present  work.  It  has  been  presented  in  the  literature  in  the  form  of  a  highly 
polished  and  formally  elegant  technique.  While  this  self  contained  form  is 
extremely  convenient  for  the  purposes  for  which  it  is  employed,  an  explicit 
exhibition  of  its  inner  mechanism  seems  to  present  certain  adventages,  and  forms 
the  viewpoint  of  the  present  report.  We  shall  first  review  the  Schwinger-  Copson 
procedure  in  a  typical  example,  and  then  develop  the  relation  to  separation  of 
variables  by  precept  and  example,  exploring  the  various  advents ges  of  the  latter 
method.  This  Involves  considerable  expository  material  which  provides  suitable 
illustrations  for  the  viewpoint  which  is  developed  in  this  report.  What  follows 
in  this  section  is  an  outline  of  the  argument  of  the  subsequent  work. 

In  Section  2  the  Schwlnger  procedure  in  its  typical  form  is  illustrated 
by  applying  it  as  an  exercise,  to  the  calculation  of  the  diffraction  by  a  straight 
edge.  This  is  essentially  the  problem  and  method  of  [9J.   In  general.  Green's 
theorem  Is  applied  to  a  contour  consisting  of  the  diffracting  screen  and  the 
circle  at  infinity.  The  knovn   Green's  function,  for  p.   simple  situation  (e.g. 
free  space)  is  employed  and  the  rave  function  to  be  sought  is  expressed  as  a  line 
integral.  When  the  boundary  conditions  are  applied,  one  obtained  by  this  means 
an  Integral  equation  for  the  discontinuity,  at  the  screen,  of  the  function  or 
(depending  on  the  boundary  condition)  its  normal  derivptive.  For  a  class  of 
problems  in  which  the  screen  or  obstacle  is  rectangular  and  Infinite  or  semi- 
Infinite,  this  integral  equation  is  of  Wiener-Hopf  type  (i.e.,  it  is  a  modification 
of  an  Integral  ecuation  first  solved  by  N.  Wiener  and  E.  Hopf  [if)  ,  (12]   The  typi- 
cal integral  equation  thus  arising  is  of  the  form 


u(x)  =  j 


,00 

g(Xo)  K  (x-Xq)  dxo  ;         x>0. 
o 
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The  kernel  K(x-Xq)  arises  from  the  Green's  function  by  evaluation  on  the  boundery. 
The  function  u(x)  is  given.  When  g(xQ)  is  found  the  solution  of  the  diffraction 
problem  is  readily  completed  by  integrations.  The  feature  of  the  integral  equa- 
tion which  is  taken  to  be  significant  is  the  difference  character  of  the  argu- 
ment in  the  kernel.  If,  In  addition,  the  range  of  integration  and  of  validity 
of  the  equation  were  (-00,  00),  the  problem  could  be  readily  treated  by  the 
application  of  Fourier  transforms,  and  their  convolution  theorem.  TMs  would 
reduce  the  problem  to  an  algebraic  one.  The  point  is  that  in  the  actu&l  situation, 
this  is  not  adequate.  However,  the  problem  caji  still  be  solved,  when  the 
customary  algebraic  considerations  are  supplemented  by  function-theoretic 
analysis  in  the  Argand  plane  of  the  (complex)  transform  variable,  the  basic 
tools  being  analytical  continuation  and  Liouvllle's  theorem.  For  convenience,  we 
refer  to  the  entire  procedure  of  section  2  as  the  Schwlnger  procedure,  which 
employs  the  Green's  function-integral  equation  fomulation,  and  achieves  a  solu- 
tion by  supplementing  "classical  transform"  methods  with  Wiener-Hopf  techniques. 

G.  Carrier*  flOj  has  shown  that  this  procedure  may  be  applied  to  other 
Integral  equations,  ^y  way  of  illustration  he  considered  a  generalized 
Tricoml  Equation 

Uyy+  y"(^-^ic2u)  =0 

and  the  boundary  value  problem  releting  to  the  diffraction  by  a  semi  infinite 
plane.  The  incident  field  was  taken  to  be  a  "plpne  wave",  i.e.  a  fimctlon  re- 
ducing to  the  latter  when  m  s  0  in  the  preceding  equation.   He  constructed  the 
Green's  function  integral  equation,  and  showed  that,  although  the  equation  is  not 
of  Fiener-Hopf  (or  difference -kernel) type,  yet,  when  a  certain  Hankel  Transform 

*  The  present  investigation  was  suggested  by  this  work,  and  by  a  private  remark 

Oi  Prof.  Carrier  that  the  transform  to  employ  in  his  integral  eouation  is  sug- 

gected  by  the  form  of  the  product  solution  of  the  partial  differential  equation 
in  ohe  form  u(x,y)  «  F(x).G(y). 
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is  applied,  (reducing  to  the  Fourier  Trensform  for  m  =  0)  an  equation  between  trans- 
forme  resulted  which  was  solved  by  a  modification  of  the  techniques  described  In 
section  2. 

The  above  work  provides  a  transition  to  the  subject  matter  of  the  present 
report.  What  we  hove  in  mind  in  this  report  is  an  exposition  of  various  relation- 
ships connecting  the  above  methods  with  the  method  of  separation  of  variables  In 
curvilinear  as  well  as  cartesian  coordinates.   In  particular,  we  present  a  character- 
ization of  the  situations  in  which  Wiener-Hopf  techniques  are  relevant  (it  must  be 
understood  that  even  where  relevant,  such  techniques  in  their  present  form  may  be 
insufficient).  This  is  seen  to  depend  on  the  form  the  boundary  conditions  take  In 
various  coordinate  systems,  a  state  of  affairs  concealed  by  the  Green's  function 
Integral  equation  formulation. 

Such  a  characterization  would  seem  to  be  useful  for  classification,  and 
integration  of  the  method  with  the  classical  procedures  employed  in  connection  with 
the  partial  differential  equations  of  mathematical  physics  with  a  view  of  the 
systematic  exploitation.   It  may  be  mentioned  (see  section  8)  that  as  s  result  of 
the  present  viewpoint  it  has  been  noticed  that  electrostatic  problems  relating 
to  an  open  finite  cone  or  conical  cup,  and  also  a  cone  punctured  at  the  vertex,  are 
now  amenable  to  treatment.  Apart  from  possible  electrostatic  applications,  such 
electrostatic  solutions  are  known  to  be  usefxil  as  the  starting  point  of  a 
variational  solution  of  electrodjrnamic  problems. 

With  reference  to  the  subject  of  integral  equations  (If  these  arise  by  a 
Green's  function  techninue  from  partial  differential  equations)  a  closer  investl- 
iretion  of  separation  of  variables  (cf,  summary  of  section  3  below)  permits  us  to 
cpecify  the  circumstances  under  which  the  state  of  affairs  is  as  exemplified  in 
[lO] ,  and  furthermore  (compare  section  6  of  the  present  report  and  also 
summary  in  the  present  section)  whether  in  addition  more  direct  and  classical 
procedures  also  apply.   In  illustration  of  these  remarks  one  may  consider  the 
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Green'E  functinn  integrel  equation  relating  to  the  conical  cup,  where  the  epplicn- 
tion  of  Mellin  transforms  together  with  Vifiener-FJopf  techniques,  would  be  suggestci 
bv  the  formulation  of  tht  boundary  condition  of  the  corresponding  partial  differen- 
tial eouotlon  problem  in  spherical  coordinates.  Section  3  conristF  of  a  review  of 
the  method  of  separation  of  Teriables,  together  with  additlonEl  articulation  of 
terminology.  The  order  of  presentation  and  the  distribution  of  emphasis  is 
here  designed  to  stress  those  aspects  which  relate  especially  to  the  pettem 
to  be  exhibited.  The  tern  "separation  of  variables"  is  reserved  for  the  re- 
presentation of  a  general  solution  of  a  linear  partial  differential  equation  by 
means  of  superposition  of  a  sum  of  particular  solutions,  the  latter  being  ob- 
tained by  writing  the  equation  in  a  certain  coordinate  system,  and  seeking 
solutions  which  are  in  the  form  of  a  product  of  a  function  of  one  of  the  coor- 
dinates multiplied  by  a  function  of  the  other.  When  a  diffraction  problem  le 
presented,  it  may  be  formuleted  as  a  boundary  value  problem  by  stating  that 
on  a  certain  boundaiy  (the  obstacle)  the  wave  function,  (or  its  normal  derivative) 
vanishes,  or  (equivalently)  that  the  scattered  part  of  the  wave  obeys  an  in- 
homogeneous  boundary  condition.   This  has  nothing  to  do  with  any  coordinate 
system,  of  course.  One  now  chooses  a  coordinate  system  In  terme  of  which  to 
represent  the  solution  by  means  of  separation  of  variables.*  In  whst  re  may 
call  the  customary  or  classical  procedure  the  coordinate  system  chostn  i  s  such 
that  the  entire  obstacle,  or.  Its  various  portions,  are  each  respectively 
obtained  by  fixing  a  particular  constant  value  for  one  of  the  coordinates ^ 
while  allowing  the  other  to  go  through  its  entire  range.  Thus  in  problems 
relating  to  a  clrcrular  cylinder,  plane  polar  coordinates  (r,e)  would  be  chosen, 
and  the  obstacle  is  given  by  r  «  constant,  while  0  <  6  <  2  tt  .  This  leads,  In 
respect  of  each  of  the  boundary  conditions  that  has  to  be  met,  to  a  simple 
boundary  value  problem.  The  problem  is  to  determine  the  coefficients  in  the 

superposition  so  ss  to  meet  certain  of  the  boundary  conditions,  and  this  is 

*  This  Is  all  said  and  done,  of  course,  subject  to  permissiveness  of  the 
sep'gr-ibility  conditions. 
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acconplished  by  che  employment  of  what  may  generally  be  called  a  '^classical 
trp-nsform  theorem".  The  particular  solutions  superposed  are  chosen  to  Beet  the 
remaining  conditions  relating  to  the  boundary  and  aire  called  proper  functions 
or  eigenfunctlons.  In  fact  the  problem  can  in  general  be  phrased  as  a 
certain  Integrcl  equation.  The  formula  for  Its  solution  is  the  trensfora 
theorem  In  cuestlon.  An  example  Is  the  Fourier  Integral  theores,  when  expressed 
as  a  pair  of  transform  equations.  From  recent  developments  It  appears  (cf.  [13j) 
that  such  a  theorem  is  uniformly  available  in  the  case  of  the  partial  dif- 
ferential equations  of  physics.  But  one  is  not  obliged  to  employ  the  suitable 
coordinate  system  in  question.  Conceive bly  the  problem  Is  formulated  in  teras 
of  coordinates  (and  elgenftmctlons)  most  suitable  for  another  problem.  For 
instence,  Cartesian  coordinates  are  best  suited  for  problems  relating  to  an 
infinite  plane,  When  the  situation  is  two-dimensional  one  may  deal  with  a  cross 
section,  e.g.  the  entire  x  axis  (y~  OJ  -oo<x<oo).  However,  one  may  be  Im- 
pelled to  employ  this  coordinate  system  for  the  case  of  a  semi-  Infinite  plane, 
y  =  0,  0<x<oo.  For  simplicity*  the  boundary  conditions  relating  to  such  an 
obstacle  may  be  taken  to  be 

u(x,y)  a  0         y  a  0        x<  0 

-|-^(x,y)  =0         y  =  0        x>0. 

Such  a  boundarj'  value  problem  1&  called  a  two  part  mixed  boundary  value 
problem.  When  the  Cartesian  eigenfunctlons  are  employed  and  superposed,  the 
determination  of  the  coefficients  Involves  the  solution  of  a  pair  of  dtial 
integral  equations.** The  use  of  the  indicated  classical  transform  theorem  is 
not  sufficient.  But  "classical  transform  methods"  when  supplemented  by 
function  theoretic  methods  or  Wlener-Hopf  technlQues  (analytical  continuation 

♦  Te   syBtematically  emit  mention  of  the  radiation  conditions  at  infinity,  and 
assume  these  to  be  met  by  the  eigenfunctlons  in  question, 

♦♦  For  the  concept  of  dual  see  cf.  [iS],  [19] . 
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and  Llouvllle's  theorem)  frequently  achieve  a  solution.  Three  part  problems  pre- 
sent additional  difficulties.  But  from  the  present  standpoint,  what  appears 
in  one  coordinate  system  as  a  three  part  problem  will  be  a  two  part  problem 
in  another,  and  a  one  part  or  simple  problem  in  a  third  system.  Such  a 
transformation  is  discussed  in  section  3,  and  illustrated  further  in  this  report. 
Most  frequently,  in  the  classical  literature,  the  choice  of  coordinates,  or  its 
implied  eq^viivalent,  the  choice  of  eigenfunctions  and  transform  theorem,  has 
been  such  as  to  correspond  to  the  formulation  of  a  simple  boundary  value  problem. 
In  these  circiunstances  there  is  a  known  duality  between  the  Green's  function- 
integral  equaxion  method  and  separation  of  variables.  One  of  the  purposes  of 
this  report  is  to  exhibit  the  persistence  of  this  parallelism,  when  the  choice  of 
coordinates,  or  eigenfunctions,  or  transform  theorem,  corresponds  to  a  mixed 
boundary  value  problem.  In  the  simple  boundary  value  problem,  the  application  of 
the  transform  theorem  immediately  determines  the  coefficients  by  an  algebraic 
eouetlon.  This  holds  not  only  for  the  eigenfunction  procedure,  but  also  for  the 
Green's  function-integral  equation  procedure  (i.e.,    the  problem  is  algebraized 
by  the  application  of  the  transform  relating  to  the  coordinate  system  and  eigen- 
functions in  terms  of  which  the  boundary  value  problem  is  simple).  The  same 
problem  might  be  treated  in  a  mixed,  two  part  manner,  i.e.  ty  use  of  eigenfunctions 
or  transform  theorem  corresponding  to  a  coordinate  system  in  which  the  problem  is 
a  two  part  mixed  boundary  value  problem.  Then  the  use  of  the  transform  alone  is  not 
enough.  This  is  the  state  of  affairs  where  Wiener-Hopf  techniques  come  Into  play. 
?.e  find  that  the  form  of  these  devices  appears  more  n'ltural  and  suggestive  in  the 
situation  arising  from  sep-^ratlon  of  variables.  But  basically,  a  complete  p-sral- 
lellsm  and  duality  still  subsists  between  the  Green's  function  and  eigenfunction 
sitiietions,  and  the  state  of  affairs  just  described  refers  to  both. 
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The  remainder  of  this  report  is  devoted  to  the  lllustretlon  of  the  fore- 
going observe tloHB  and  accordingly  to  the  carrying  over  of  the  Wlener-Hopf 
techniques  Into  the  eigenfunction  viewpoint  where  appropriate.  In  section  U 
we  carry  out  the  eigenfunction  treatment  for  the  semi-infinite  plane  (i.e. 
straight  edge)  problem,  when  viewed  as  a  two  part  problem.  This  is  the  couter- 
part  cf  the  Schwinger-Copson  treatment  of  section  ?.  The  greater  freedom  of 
viewpoint  reveals  some  features  reletlng  to  uniqueness,  which  ^re  here  discussed. 

In  section  5  a  more  complex  problem  of  the  same  type  is  discussed.  The 
reflection  or  diffraction  properties  of  an  infinite  staggered  array  of  semi- 
infinite  plates  have  been  treated  by  Carlson  &  Heine  £nj  by  whet  may  be  des- 
cribed as  an  orthodox  application  of  Bchwinger's  methods.  Since  the  geometry 
here  involved  seems  far  from  separable  it  seems  of  interest  to  show  how  the 
present  viewpoint  applies  there  too,  and  hence  section  5  is  devoted  to  a  treat- 
ment of  this  problem  by  separation  of  variables. 

The  above  sections  relate  to  the  employment  of  Cartesian  coordinates, 
and  Involve  the  standard  Fourier  transform.  In  the  subsequent  sections  polar 
coordlnrites  are  in  cuestion.  Thus  in  section  6  the  problem  of  the  semi-infinite 
plane  is  examined  from  the  standpoint  of  plane  polar  coordlnntes.  In  terms  of 
these  coordinptes  the  boundary  value  problem  becomes  a  one  part  problem.  The   eigen- 
functions  are  Ber-sel  functions  J  (kr)e-      and  the  transform  theorem  dealing 
with  the  superposition  of  these,  veryine  the  index  or  order  V  ,   has  been  devised 
by  N.  N.  Lebedev  and  employed  by  Lebedev  and  Kontorowitsch  ^5]  in  a  separation 
of  variebles  (or  eigenfunction)  treatment  of  the  problem  in  a  classical  manner. 
Here  re  show  how  polnr  coordin^.tes  a-^e  employeble  in  a  parallel  manner  in  the 
Green's  function  -  integral  eouatlon  formuletlon  of  section  2  in  that  the  appli- 
cation of  the  Bessel  transform  reduces  the  integral  ecuatlon  to  an  algebraic  one. 
The  integral  eruatlcn  of  section  2  Is  thus  a  true  convolution  ee  far  as  the  Bessel 
transform  is  concerned,  and  Wiener-Hopf  techniaues  are  not  necessarj-. 
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This  illuntrateB  the  possibility  also  that  one  and  the  same  Integral 
eouatlon  may  be  Bolved  by  Tfiener-Hopf  Techniques,  and  In  addition  by  more 
direct  methods  employing  a  different  transform.  This  reflects  a  change  to  a 
more  natural  coordinste  system  in  the  partial  differential  equation,  a  change 
to  which  the  integral  equation  is  insensitive. 

In  section  7  a  related  problem  Is  considered.  We  consider  a  static 
problem  relating  to  an  infinite  ribbon  whose  cross  section  occupies  the  x  axis 
from  X  =  0  to  X  =  1.  Laplace's  equation  in  two  dimensions  is  the  equation  we 
deal  with.  A  formulation  in  the  manner  of  section  2  indeed  leads  to  an  integral 
eouatlon  with  a  difference  kernel,  but  with  limits  of  integretlon  (0,1).  Direct 
further  imitation  of  Schwlnger's  Fourier  transform  procedure  fails.  But  we  may 
view  this  situation  as  arising  from  the  fact  thpt  use  of  the  Fourier  trs.nsform 
corresponds  to  the  employment  of  Cartesian  coordinates,  in  terms  of  which  we  have 
a  three  part  problem.  The  transformation  to  polar  coordinates  leads  to  a  two 
part  problem;  however,  the  radial  eigenfunctions  are  powers  of  r,  and  the  Mellin 
transform  Is  suggested.  Te   thus  disregard  the  difference  kernel  aspect  of  the 
matter.  A  separction  of  variables  in  polar  coordinates  leads  to  a  solution, 
when  the  Mellin  theory  is  employed,  in  conjunction  with  Wiener-Hopf  techniques. 
It  is  clear  that  the  Green' £  function  -  integral  equation  may  also  be  so  treated. 

In  section  8  the  problem  of  the  charge  distribution  on  a  disc  Is  brought 
within  the  scope  of  methods  of  this  report.   This  problem  has  been  treated  by 
Titchmarsh  j_l6J  from  the  standpoint  of  cylindrical  coordinates.    In  that  treat- 
m<=rit,  8  3et  of  dual  Jntepral  ecuotlons  wrlse  from  use  of   the  cylindrical  coordl- 
nates.   These  ecuatlons  are  trarfformed  so  as  to  f?lve  rise  to  a  Mellin  trnnsform 
problem  which  is  solved  infeniously.   The  present  treatment  eeeir.s  to  clarify  the 

P..   E.  Heins  has  recently  reworked  the  solution  of  this  pair  of  dual  integral 
eniatlons  using  f/ellin  trsnrifonnE  and  Fiener-Hopf  Technioues.   cf-^igj. 
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phenomena.  The  problem  Is  treated  In  spherical  polar  coordinates,  as  a  tw>  p&rt 
problem  of  a  type  uniform  with  that  of  the  preceding  section.  The  transition  fron 
cylindrical  to  spherical  coordinates  is  in  fact  fornally  a  transforastlon  fro« 
"cartesians"  to  polars.  This  use  of  spherical  coordinates  leads  to  the  eigen- 
functlon  r~  P^.  •,  (±  cos  9)  and  this  directly  to  Mellin  trensforn  theory.  It  is 
furthermore  readily  generalizable  to  problems  relating  to  a  truncated  or  finite 
cone.  A  solution  to  the  latter  has  been  obtained.  It  is  touched  on  here  and 
constitutes  the  subject  of  a  separ«.te  report. 
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2.  Illustration  of  the  Green's  Function  Integrnl  Fcu?tlon  Method.  Dlffr'^.ctlon 
by  a  Straight  Fdge. 

In  the  present  section  we  wish  to  illustrate  the  method  of  Schwinj^er 
(and  Copson)  by  the  treetment  of  diffraction  of  «i  plnne  wave  by  n  strpl^ht 
edsre,  or  peml-lnfinite  plane.  The  application  of  Green's  function  snd 
Green's  theorem  leads  to  an  integral  eountion  of  modified  Wiener-Hopf  type, 
which  is  solved  by  function-theoretic  methods.  This  is  an  extremely  simple 
example  of  the  class  of  problems  solved  by  Schwinger  in  this  manner,  and  is  given 
here  in  order  to  present  the  essential  features  of  the  method  as  outlined  in  the 
introduction. 

The  edge  of  the  semi -infinite  plane  screen  is  taken  as  the  z  axis  of 
a  cartesian  coordinate  system.  The  screen  occupies  the  right  half  of  the  x,z 
plane,  i.e.,  y  »  0,  x^O.Csee  flgurellWe  will  assume  th?t  the  ""irectlon  of 
motion  of  the  incident  wave  lies  in  a  plane  perpendicular  to  the  z  axis. 
Then,  by  symmetry,  the  problem  is  a  two-dimensional  problem,  i.e.,  conditions 
do  not  depend  on  the  value  of  z.  Te   will  use  acoustic  terminology.  Let  u 
represent  the  velocity  potential,  p.nd   let  u  represent  its  sp  ce  dependence. 
Then  u  =  n  (x,y)  e~^*^  with^Os  kc,  k  »  -jr".  Here  A  is  the  free  spsce  \«Tve 
length,  and  c  l-^  the  velocity  of  sound.  We  consider  *?  plane  wnve  incident 
upon  the  screen,  and  write 


\ 


(?.l)  u  =  u  + 


where 


_     ik(x  cos  9g-«-  y  sin  Q^) 
o  ~ 
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represente  the  Incident  wave,  snd  where  tL  represents  the  scattered  irave.  The 
geometry  is  indicated  in  the  adjacent  flgxire. 


Direction    of 
incident    wove 


screen 


■^  X 


Fig.  I 


The  boundary  value  problem  may  be  formulated  as  follows: 


(2.3) 


u      4-  %y-t-   k  u  »  0 


k  = 


=  k,-»-ik2;     0<k2<<l 


(2.-C) 


_e  u 
ay 


=0;     X  >  0 


y  -  0 


(2.5) 


U-   — ♦     0  at  oo,  away  from  the  screen. 


(2.6) 


u  regular,  and,  in  partlculer,  continuous,  in  the  exterior 

of  the  screen. 


It  will  be  noted  that  a  small  pjositlve  imnginary  part  has  been 
assigned  to  k.   This  corresponds  to  e  slight  degree  of  absorption  in  the 
medium,  and  is  extremely  convenient  mathemotlcally. 
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One  now  Introducee  the  free  space  Green's  function.  TMb  ie   defined 
as  a  solution  C(x,y;xQ»y  ),  possessing  a  logarithmic  singularity  at  the  point 
(xQ,yQ),  and  vanishing  et  Infinity  .  The  possession  of  a  logarithmic  singu- 
larity ^-ives  it  the  cherecter  of  an  oscillating  line  mass,  in  analogy  with 
the  function  log/  (x-x^)  +.  (y-y  )^   which  represents  a  line  mass  at  (xQ,y  ) 
In  two-dimensional  potential  theory.  Such  a  line  mass  may  be  regarded  as  a 
density  distribution,  Tdth  a  density  which  is  infinite  at  (x  ,y  )  and  zero 
elsewhere.  The  wave  analogue  of  Poisson's  equation  for  a  density  distribuoion 
is 

^^  u  +•  kru  s  p 

where  p   is  the  density.  The  analogue  of  the  total  bbes  ie  \\P   dA  where  the 
integral  is  over  the  entire  plane.  This  must  be  equal  to  1.  Hence  f  is 
zero  everynhere  except  at  (xQ,y  ),  where  it  ie  infinite,  and  its  integral  is  uni- 
ty. Such  a  function  it  denoted  ty  the  symbol  -  5  (^-^)  5  (y-yo)»  and  called 
the  Dirac  delta  function.  Thus  G  obeys 


^^a(x.y;x  ,y  )-^  k^G 


O'  0 


-  6  (x-x^)  6  (y-y^) 


The  solution  is  known  to  be 


H 


(1) 


G(x,y:Xo,yo)  =  -;^  "o^^^  [k  /(x-x^)^  +  (y-y^)^ 
is  the  Hankel  function  of  the  first  kind. 
Now  Green' £1  second  identity  is 


rr(uAY-vAu)  dA  =  i    (u-^-v.aJl) 
•'•^A  Jc  «3n     an 


dJS 


,  where 


where  C  ie  a  closed  curve  bounding  a  region  A,  and  where  ^ —  denotee  the 


*  In  consequence  of  the  asstimed  imaglnnrj"  part  of  k,  it  is  enough  to  require 
vanishing  at  oo . 
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norre-.l  derivrtlve.   If  u  is  our  potential  function,  and  v  la  6(T,y;x  .y  ),  the 
Green's  function,  and  If  C  Is  the  Indented  circle  of  large  radius  In  the 
appended  figure. 


screen 


C  =  Ci-hC2 


then,  using  the  equation 


Fig.  2 


uaG-Gau  =  u  [  -k^G  -    S  i^'\)   6  (y-y^j)  ]     -G(-Ic^u)  =  -  u^(x-x^)  5(y-y^) 


we  have 


u(x,y)  =      J      tG(x,7;xo,yo)  -|-^  (xo»yo)  "  u(«o»yo)  "g^  (5c,yj3Co,yo)]  ^S,  . 


If  we  now  let  the  radius  of  the  circle  go  to  Infinity,  the  contribution 
froB  Cp,  as  a  consequence  of  the  plane  wave  behavior  of  u,  is  u  .  We  now  let  C^ 
approach  the  screen  and  use  the  boundary  condition  -^■-  =  0  on  the  screen,  to 
obtain 


(2.7) 


TOO 
u(x,y)  = 


i^-|r  (''yj'^o'yo) 


ay 


dx  +  e^^^'^  °°®  ®o*  y  ^*"  ®o^ 


o 

y  =  o 


where  [u]  denotes  the  difference  in  the  ?alue  of  u  on  the  upper  and  lower  sides 
of  the  screen. 

In  the  integrand,  consider  the  factor  -^ —  .  Since  G  is  logarithmic,  its 

ay© 

normal  derivative  represents  a  dlpole.  The  integral  represents  so  to  sp>eak  the 
potential  of  a  double  layer,  whose  density  consequently  is  proportional  to[u(x  )J, 
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ThuB[u(x  )j  repreBents  the  die  continuity  of  u  on  pessage  throii^h  the  l:;yor.  We 
may  Identify  the  Integral  ulth  u,,  the  scattered  wave.   If  »e  carry  out  the  dif- 
ferentiation we  find  that  u   containB  a  factor  y,  and  Is  an  odd  function  of  y. 
The  discontinuity  of  u  Is  also  that  of  u^ ,  since  u  is  continuous.  On  the  other 
hand,  u,  is  continuous  at  y  =  0,  for  x  negative,  since  u  and  u  are;  and  hence '^ 
u,  should  vanish  on  the  negrtive  x  axis.  This  Is  seen  to  be  the  case  since  y 
i£  a  factor  of  u, . 

If  we  now  apply  the  boundary  condition  .P.  ^     =  Oj  x  >0, 

9  7J 
and  denote  [u]  by  f (x  ) ,  we  have  X  "  o 

(2.8)    0  =  ik  sin  e^  e^'^'^  °°^  ®o  -V  |    f(x^)  ^^^  (x,0;x^,0)  dx^  j  x  >  0. 


Since  the  argument  of  G  is  /  (x-x^)^  -4-  (y-y^)   ,  one  finds  thnt  the  argument 

2 
of  ^  G     (x,0;x-,0)  is  |  x-x  I,  where  the  absolute  value  is  taken  by  virtue  of 

^ydyo     °       ° 

the  feet  that  the  argvunent  represents  a  distance,  and  hence  must  be  positive. 


Thus  our  integral  equation  is  of  the  form, 


,.r 


(2.9)  0  =  h(x)+  I   Hx^)   K  (U-XoO  d^o  »  *  ^°' 

o 


where  K  replaces 


9^G    -   K)J^Miii2^) 


Attention  Is  no*  fixed  on  the  argument  of  the  kernel  K.  Since  this  Is 
j  x-x  I  we  are  reminded  of  the  convolution  theorem  for  Fourier  integrals.  Es- 
sentially this  states  that  if 
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-oo 


•/ 


-oo 


oo 


K(xJ  e"*"^"^  dx 


-oo 


are  the  trenefomiB  of  f(t)  and  K(t)  respectlvelj,  then  the  transfom  of  the 

/-  oo  _    _ 

expreeelonl    f(x^)  KCx-x^)  dx^  ie  f (« )  K(ar).  Hence  if  our  equation  were 

J  -co 
valid  for  all  x,  and  the  renge  of  integration  were  fron  -oo  to  -»•  oo ,  we  could 

te.ke  Fourier  ti^nsforms  and  obtain 

0  =  h(cy  )  +  f((x).  K(o^)   . 

Then  f(cx  )  =  -  ^L2Ll  ,  and,  ly  the  Fourier  trensforu  theorem 
K((X) 

°    ^^  j  ^TT  J    K(oc) 

-oo  _oo 

Thus  fCxp)  would  be  readily  obtained.  The  special  advance  of  the  method 
under  consideration  is  that  the  equation  is  solved  despite  the  fact  that  the 
range  of  integration  is  o  to  oo  and  the  Fourier  transform  theorem  carjiot  be 
applied  directly.  One  extends  (2.9)  to  read 


0(x)  =  h(x)+  I 


oo 
(2.10)        0(x)  =  h(x)+  I    f(x„)  K(»x-Xol)  dx„  ;  -oo<x<cao 

-oo 
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wlth  the  convent ions  that 

0(x)  =  0:  X  >  0 

(2.11)  I   h(x)  =  0;  x<  0 

This  is  equlvPlent  to  saying  that  while  according  to  (2,9)  the  ■value  of 
the  integral  involved  is  knourn  as  -h(x),  for  x  >  0,  its  unknown  values  for  x 
negative  are  given  the  name  0(x) .  The  Fourier  transform  can  now  be  applied  and 
we  have 

(2.12)  ^(«)  =  h(£x')-i-  ?(«)  K(CX) 

but,  now  in  this  equation  ire  have  two  unknown  ftmctions,  0(«)  and  f  (  oc) .  How- 
ever both  can  be  deduced  slnsultaneously. 

Actually,  according  to  Fourier  integral  theory,  one  has  to  Investigate 
doraainE  of  regularity  of  0,  h,  f ,  and  K  as  functions  of  o*  ,  where  the  latter  is 
regarded  as  a  coir|.lex  variable.  The  equation  (2.12)  only  holds  in  a  common  domain 
of  regularity  if  such  exists.  We  make  use  of  the  complex  character  of  k  in  this 
investigation. 

In  cur  case  h(i)  =  0  for  x  <  0,  and  h(x)  =  ik  sin  9^  •  e^^  '^°^   ®o  for 

X  >  0.  Thus  i(cx)  =  ik  Bin  %    (^   e^^^  ^°^   ®o  '  "^  ^dx.  The  integral  con- 

J  o 

verges  imiformly  when 

Im(cx)  <  Im(k  cos  60)  ^1  ^2  cos  Qq       , 

where  kj  Is  the  ImaglnRry  part  of  k.  Hence  h(oi)  is  regular  in  a  lower  half 
plane  which  includcp  the  real  axis,  and  is  there  represented  by 
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_        ksin  e^ 
(2.13)  h(<x)  = 


or-  k  COB  6g 


Row  f  (x^)  is  proportional  to  the  Talue  of  ^n)  on  the  obetacle,  and  one 

i  K3Cn 

assunes  on  phyaieal  grounds  that  f  (x  )  '^9  9   where  M  has  a  positlTe,  is- 

aglnar7  part.  Then 


00 
f(«)  5  f   f(Xo)  b'^*^^  dx^ 


is  regular  in  a  lower  half  plane,  including  the  real  axis,  where 

lm(o<.)   <   iMfH) 

Using  the  asymptotic  form  of  the  Hankel  function  type  of  kernel,  i.e., 

H^  '  "^  e       for  X  large  ^.nd  negative,  and  using 
o 

equation  (2,10),  we  find  that  i^(x)  is  asymptotic  to  e~     for  x  large  and 
negative,  where  k  is  a  complex  number  with  a  positive  inaginary  part.  Using 
this  fact,  and  using  the  definition 


?(  «  )  =1    0(x)  e"  **  dx 


«)=r    0(x)  e' 
-00 


we  find  that  ^(o<  )  is  regulnr  in  the  upper  half  plane,  Tm(o«  )>   -  ImCk  ), 
Also,  it  is  known  that  K(«  )s   -  7k^-o<,  regulnr  in  the  strip 

|lm^)|<  ImCk). 

Thus,  surrounding  the  real  axis  is  a  common  strip  of  regularity  of  all 

the  transforms,  and  in  this  strip  equation  (2,12)  hold3,  and  in  our  case  takes 

the  form 
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One  notes  that  the  function  ^  i^  -'  of      y   which  arises  from  the  Green's 
function  kernel  of  the  Integral  equation,  can  be  factored  into  two  functions 
y  k  -  o<   and^k  t «  ,  which  are  regular  in  upper  and  lower  half  planes  respectively. 
This  is  quite  generally  the  case. 

Tfe  use  the  subscripts  -♦•  and  -  in  what  follows,  to  indicate  regularity 
in  upper  and  lower  half  planes  respectively.  From  (2.1A)  we  have. 


XL^  «  f-jL i    JL^}lIo_  ^    if  (PC  )  n^ 

x/kV^"    yirr^J    o^ - ^  cos eo         2     ^^' 


ex. 


or 


Hgg ,  .^'^° .^  f -.^A^ .      1      \^  " '^° "o /     i„vL>Ll 


\y/k7^      /ktkcose^   )   cx-kcoseoV^^os  9 J   2    '^- °^ 


The   first  term  on  the  right  no  longer  has  a  singularity  at  a  s  k  cos  Qq,   but  is 
regular  in  the  upper  half  plane.  Hence,  if  we  write. 


^^,-,  ^-^ A,,.J\.  .-■■■.  ■■  ,„.  q ». 

/k-l-o<   -/ktk  co88o/~  -v^kfk  coseJ((o<  -kcos9Q) 


then  in  this  equation  the  left  side  is  regular  in  an  upper  half  plane,  while  the 
right  side  is  regular  in  a  lower  half  plane.  Since  these  half  planes  overlap, 
the  equation  (2.15)  may  be  regarded  as  defining  a  certain  entire  function.    The 
left  and  right  hand  sides  give  its  representations  In  upper  and  lower  half  planes. 
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The  final  device  in  th5s  procedure  is  the  determination  of  the  entire 
function  in  question.  This  will  be  accomplished  by  a  fora  of  Llouvlll^e 
theorea,  which  states  that  If  an  entire  function  of  of.    is  known  to  behav.* 
like  o<   as  o<  — f  oo.  It  mist  bs  a  polTnomlal  of  degree  n.  In  the  present 
ease  one  can  study  the  growth  of  the  entire  function  in  the  upper  and  lower 
half  planes  by  studying  its  respective  representations.  One  finds  that  the 
function  behaves  like  ex  P  where  p  is  negative.  Hence  It  must  be  at  most  a 
polynomial  of  negative  degree.  It  must  therefore  be  identically  zero. 

To  study  the  growth  one  investigates  the  beV^vlcr  of  the  transforms 
involved.  We  have,  in  the  lower  half  plane,  o<  -    "^x"*  ^  ^2*   *•***  ^^  2 
negative. 

f(c.)|  i=.    I  °°  |f(xo)  e"''''^!  dx  -  r      If(xj)   e-  I' ^l^°  dx^ 


[C)\    ^  |f(xo)  e"   '^l  dx^£      jf(x^)j 

Jo  '      •'  o 

In  order  to  have  been  able  to  write  (2. 7)f (Xq)  must  be  integrable  at  the 
origin,  and  thus  fCxg)  may  be  assumed  to  behave  near  x  =  0,  like  x^  , 
where  p  >  0. 

Now  it  is  clear  that  the  largeness  of  o*   produces  a  strong  decay,  so 
th5.t  the  integral  vanishes  as  «  — >  oo  with  la  {.o< ^   negative.  But  further- 
more the  effect  of  f (x^)  is  negligible  almost  everywhere,  due  to  tAis  fact,  with 

-i  o<  X- 
the  exception  of  the  neighborhood  of  x^  «  0,  where  e     -o  ^^  j^  inde- 
pendently of  «  .  Hence,  the  behavior  of  f(o<)  can  be  well  approximated  by 
using  the  form  of  f (x^)  at  the  origin  .  Thus 


-;:?i- 


Dsing  the  substitution      |oc  I  x^  *  8,  we  obtain 

f  (cK  )  '^  (X         t   as  ex  — >  CO;   Iffl(or)<  0  • 
Similar  analysis  leads  to 

^^  («)  ~  ©<■  ~^,      erf  -^  oo;     Im(oO  ?  0, 

if  we  assioffie  0(x)  ■^  (-x)     as  x  — >   -  0.  Again  too,  q  >  0. 

Hence  the  asymptotic  form  of  the  right  side  of  (2.15)  is  cC     , 
or  like  (x  ~  »  whichever  is  greater.  For  the  left  side  we  find  the  be- 
havior ex  ~^~*   .  Since  p  and  q  are  positive,  our  entire  function  is  a 
polynomial  of  degree  less  than  ^,  and  hence  It  Is  of  degree  zero,  i.e.^a 
constant.  That  this  constant  must  be  zero  may  be  concluded  from  the  fora 
of  c<  ~^~'^  since  q  is  positive. 

Hence  from  (2.15),  since  the  left  and  right  sides  must  each  be  zero, 
we  infer  that 

2ik  sin  0^  , 


(2.16)  f(«) 


V  k*k  cos  e^  Vr^     '        o^  -  Ic  cos  9q 

Hence  in  fact     f(o<)    r>^  c<  '^         so  that  p  =  -»■  3/2.     Thus 

(2.17)  f(x  )~  X  asx^— >-l-0 

^         '  o  o  o 


Also 


oo 


-OO 


f(V  =    -r^     I  f(o<)   e^     '^  d« 


or 

/OO 
e^°*'°     dcx 
_^      (o(   -k  co8  9o)/-irr^ 

Since  fi^^^)   is  [u],   this  can  be  in  principle  inserted  into   (p.y)   to  pro- 
vide an  explicit   if  complicated  ejcpresslon  for  u(x  y). 
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3.  The  Method  of  Separation  of  Variables. 

In  this  section  ne  vlsh  to  review  and  dlseuciB  the  aethod  of  separation 
of  varieblea  in  a  general  manner,  with  a  view  to  characterising  explicitly 
the  domain  of  applicability  of  the  classical  methods,  and  the  Wlener-Hopf  Tech- 
niques, Fe  shall  also  introduce  soae  dictinctions  in  terminology. 

Consider  the  typical  linear  second  order  partial  differential 
equation  in  two  variables  which  arises  in  physical  problems.  This  is 
customarily  an  invariant  differential  equation,  i.e.,  the  ouantitles  Involved 
have  intrinsic  geometrical  or  physical  significance  apart  from  the  coordinate 
system  in  which  they  are  expressed.  To  make  the  ideas  concrete  we  may  keep 
the  equation  a  U4-k^u  =  0  in  mind,  since  the  applications  are  chiefly  to 
this  equation.  Suppose  we  realise  the  differential  invariants  in  terms  of  a 
certain  coordinate  system  ^  ,  1  ,  where  ^  ~  ^  (^c,/) »  '^  =  1  (3r,y),  and 
where  x,y  are  cartesian  coordinates.  The  differential  equation  if  of  elliptic 
type,  may  then  take  the  form 

(3.1)  u^^  ^   u^^  +   F(^  ,^]  )u  =  0 

Now  suppose  that  there  are  solutions  of  the  form  u  =  ^(  T  )  y^{  ?l  ). 
If  this  is  the  case  the  differential  equation  is  said  to  be  separable  in  the 
coordimite  system  f  »  ^»  "•  'then  obtain,  after  substituting  in  (3.1)  and 


♦  The  case  of  polar  coordinates  does  not  at  first  appear  to  be  of  this  na- 
ture, for  inLaplace*s  equation  we  have  ^  -"•■  "Z  U-"*"  "j?  ^na   *  ^*  ^* 
if  we  introduce^  =  log  r  we  have  u..  +•   *?  qq  "  ^*  ^°  general  we  ob- 
tain the  canonical  form  in  question  by  choosing  ^  »  *l   os  conjugate 
harmonic  functions  of  x  and  y,  i.e.,  ?  +  i f^  =  fix  +  ly)  where  f  is  an 
analytic  function. 
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dl riding  by  u  =  f^*  y' 


(3.2) 


0  if  )      -¥'(;7)*'^M'^>-° 


?  )        y 

without  giving  a  complete  discussion  of  possibilitleB,  we  may  say 

that  (3.1)  separates  if  F(  5"  »  *?  )  =  Fj^(  f  )  4  ^2^ 'X  ) »  ^^^^^  ^i   *n^  ^2 
nay  be  constants  or  zero.  One  finds 


g"(  f  )  +  F  (  f  )1  +  [  fin.)    +  F  (  A)  ) 

Ttf)  *  ^M  ^  ^         ftT)   *  ^2^  'I  ^ 


=  0 


Since  the  first  bi^cket  depends  on  ^    alone,  and  the  second  on  r^    alone  the 
Inference  is  that  each  must  be  constant,  and  the  sum  of  the  two  constants 
must  be  zero.  Ve   may  thus  achieve  a  solution  u(^,'7,/l)=  0  (  T)  1^(  *?) 
if  0,  J^are  so  chosen  that 


(3.3) 


These  ordinary  differential  equations  are  known  as  differential  equations  of 
the  Sturm  Llouvllle  type. 

Since  the  constant  A  ^^7   ^  chosen  arbitrarily,  and  since  the  dif- 
ferential equation  is  linear,  a  sum  of  terms  u(  ^  ,  »2  >  ^  )   corresponding 
to  various  values  of  A  is  also  a  solution,  and  we  may  thus  expect  to  build 
up  a  general  solution  by  setting 


(3.^) 


u  =    j a( A  )  0^  ( o  f,in)  ^^ 
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We  mny  view  the  integral  as  a  cent our  integral  In  the  complex  A  plane. 
The  contour  is  as  yet  unspecified.  If  the  contour  is  closed  and  contains 
poles  of  the  integrand,  say  A    A  .......  A _,  and  the  integrand  is  other- 
wise regular,  we  expand  by  the  method  of  residues,  and  obtain  a  series 

^    ^^I^X       ^  O  ]^   (  '7)«  '"^en  the  contour  is  "closed"  in  the  infinite 

n  n 

part  of  the  plane  the  number  of  poles  may  be  infinite,  and  we  hrtve  a  series 

representation,  which  is  thus  included  in  the  integral  representation.  Hhen 

the  singularities  of  the  integrand  are  of  a  more  complex  type,  such  as  branch 

points,  we  may  perhaps  view  the  branch  line  as  a  continuous  distribution  of 

poles,  and  thus  obtain  an  integral  representation  instead  of  the  -repreEentation 

in  series.  To  fix  the  ideas,  we  may  point  out  that  the  series  corresponds  to 

the  well  known  Fourier  series,  and  the  integral  to  the  Fourier  integral.  This 

method  of  representation  of  a  solution  of  a  p&rtial  differential  equation  as 

a  sum  of  special  solutions  i^»  (  C  )  ^|  (  *7  )  *®  will  call  as  is  customary 

the  method  Separation  of  Variables. 

>    All  the  foregoing  has  no  reference  to  the  subject  of  boundary  value 

problems.  Usually  however,  we  are  asked  to  find  in  a  certa  in  region  a  solution 

which  takes  on  certain  values  as  we  approach  the  boundaries  of  the  region. 

Frequently,  the  method  of  separation  of  variables  is  employed  in  such  cases. 

The  Invariant  differential  equation  is  written  in  a  suitable  set  of  coordinates 

f  ,  >|  ,  What  is  meant  by  suitable  is  the  following:  one  chooses  (if  possible) 

C  »  *7  »  to  be  such  that  the  equation  separates,  and  such  that  boundary  is 

composed  of  parts  along  which  ^  and/or  \\  are  constant.  The  equations 

%  (x»y)  -  constant 
(3.5) 

^  (3f»y)  -  constant 

represent  a  network  or  family  of  curves  in  the  x,y  plane.   If  the  boundary  of 
the  region  in  question  is  included  in  this  netuork,  part  of  our  aim  is  achieved. 
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^^ie=b,£.Q 


Fig.  3 


>■,  -r~  \   may  be  prescribed.  In  the  first  case  fie   have,  for  example, 


The  boundary  condition  mey  be  given  as  follows:  as    s  ■ — *  ?©  » 
a  constant,  u(  ^  ,  *7  )  —»  ^i(  C  i  *^  )  =  f(  *^  )»  where  f (  v^  )  is  a  given  function. 
Alternatively, 
the  equation 
(3.6)    u(  5  ^,  n  )  =  f(l  )  =   (a(  A)  !^^(  ^^)  3^(  *?)  dA   ^i  -^  ^  ^  ^2  , 

where  we  suppose  that  t^  varies  from  ^t  to  V^2  ^long  the  curve  f  -    ^    .     In  the  second 
case  we  have 

^^•"^^   ^|=e(^2)  =  I  a(A)  ^^    Y'^iYi)  dA,^^  ^  yy,    ^  nz     > 

If,  in  (3.6)  f (  >^  )  is  given  over  the  range  V|^  .<i.  ^^  <■    ^    ^,   the  equation  is  an  in- 
tegral equation  for  a(  A  )  •  Ii^  many  cases,  the  explicit  solution  of  this  equation  is 
well  known.  In  other  cases  it  can  be  obtained  by  certain  recently  devised  procedures. 
Our  purpose  is  not  the  discussion  of  these,  and  we  refer  the  reader  to  Titchmarsh  [  13^ 
for  a  me.themstlcal  treatment  of  the  subject.  T^e  merely  wish  to  recall  that  in  many 
cases  the  solution  of  equations  like  (3.6)  or  (3.7)  consists  of  a  simple  integral  ex- 
pression for   8(  A  )  .  For  example,  in  the  case  of  equation  (3.6)  we  might  have, 

(3.8)  ^^  ^  ^  f^A^  ?o^  "   r  f  (  r.  )  ^*  (  n  )  L  (  K?  )  d  Kj 

/*  ■'A 

where  ^A*  is  a  linearly  Independent  solution  of  the  same  equation  as  that  for 

1^  (l^),  and  LCtrj)  is  a  prescribed  -^ei^ht  function.  Let  us  set  a(  A.  )0(  ^  )=A(/V). 
Then  f (  >o  )  is  called  the  V  transform  of  A(  A  ),  and  A( A  )  is  called  the  inverse 
transform  of  f (  >^  ).  The  pair  of  formulae  constitute  a  transform  theorem.  A  well 
known  case  is 
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(3.9) 


VTjf  J 


-oo 


-oo 


-oo 


the  theorem  of  FoTirier  transforms.     Of  course  the aq>res8ion  for  A(  y^  )  way  In- 
volve a  simple  pole  integrand,  and  lead  to  a  pair  of  formulae 

f(v,)=  II A_  fjn) 

(3.10) 


A„=      jfin)  fj^   )  L  (ri)  dri 


In  this  familiar  case,  the  second  of  equations  (3.10)  is  the  formula  for  the  »-th 
Fourier  coefficient  of  f (  Kj  ),  and  the  range  of  integration  is  the  interval  of 
orthogonality  of  the  set  of  functions  y^. 

Suppose  now,  however,  that  instead  of  either  (3.6)  or  (3.7)  we  have  the 
following  set  of  equations 


u(^  ^,»x)  =  f(n)  =  /*(  ^)  <^;^(  ^o^  )^^<»?)  AA  »  n^^n^  n 


(3.11) 


d 


f(^o>^)  =  g(n>  =  j»(^  )i^x(l^  ]^(ri)dn,  r{  ^h^^: 
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ThuB  a(  A  )  obeyp  a  pair  of  d\ial  Integral  equations.  I.e.  a  dlffei-ent  equation  in 
each  part  of  the  range  ^,  ■^  1  '^  1o«  Equation  (3.8)  cannot  be  resorted  to  be  a 
solution.  In  other  T»ords  the  known  transform  theorem  is  not  p.ppllcable  he-^e.  In 
the  discrete  ca?e,  T-here  we  are  solving  for  a  Fourier  coefficient  this  msy  be 
simply  attributed  to  the  fact  that  the  range  of  validity  of  either  of  equntions 
(3.11)  ir.  not  the  range  of  orthogonality  of  the  ftmctions  ^(  to).     We  note  that 
the  present  situation  arises  when  on  pert  of  the  boundary  f  =  ?«  (conetant), 
the  function  is  given,  while  on  the  other  part  the  normal  derivative  is  given. 
This  is  called  a  Mixed  Boundary  Value  Problem.  On  the  other  hand,  the  more  for- 
tunate situation  previously  discussed  is  one  in  which  the  boundary  ^  =  T   is 
not  divided  into  several  parts,  but  rather  the  function  u  or  its  derivative  alone, 
are  specified  all  along  the  curve.  This  we  will  call  the  Simple  Boundary  Value 
Problem.  We  may  also  envisage  the  case  where  the  coordinate  line  f  =  f o  ^^ 
divided  into  3  or  more  parts,  various  combinations  of  u  and  -—-^   being  given  on 
each. 

It  Is  clear  that  the  solution  of  a  simple  boundary  value  problem,  in  the 
case  of  a  single  integral  equation,  is  equivalent  to  the  possession  of  a  transform 
theorem  and  is  directly  solved  where  this  is  known,  as  it  is  in  most  cases.  And 
the  range  of  known  cases  can  be  readily  extended  using  for  instance  the  methods 
explained  in  Titchmarsh  [l3]  .  In  the  present  report  we  wish  to  emphasize,  rather, 
the  case  of  a  two  part  boundary  value  problem.  We  wish  to  call  attention  to  the 
fact  that  it  furnishes  the  circumstances  in  whi^h  T'iener-Hopf  techniques  come  into 
play.   From  the  present  standpoint  we  can  say:  when  a  problem  has  been  formulated  li 
a  separable  coordinate  system  in  terms  of  which  it  is  a  one  pnrt  proMi^'m,  the  cor- 
responding transform  theorem  suffices  for  its  solution,  both  in  the  separation  of  va- 
riebles  or  in  the  Green't  function  integral  equation  formulation.  When,  on  the  other 
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hand,  separation  in  the  coordinates  system  In  question  leads  to  a  two  part 

problem,  additional  function-*theoretic  techniques  are  called  into  play,  and 

may  succeed  in  furnishing  a  solution.  These  apply  to  both  the  Green's  function 

Integral  equation  and  to  the  method  of  separation  of  variables. 

The  three  part  boundary  value  problem  does  not  appear  directly  amenable 

to  treatment,  but  the  previous  considerations  suggest  the  following  device. 

Suppose  we  choose  a  certain  coordinate  system  leading  to  a  three  pert  boundary 

value  problem.  We  may  then  change  to  a  new  coordinate  system,  or,  what  Is 

equivalent,  introduce  a  transformetion  of  coordinates,  such  that  the  boundary 

curve  Is  included  in  the  coordinate  mesh  but  is  characterized  by  two  different 

values  of  f  ,  say  ^  =  E,    ^  "  ^  1»  ^^   ""^y  thus  achieve  the  "Wlener-Hopf" 

state  of  affairs.  The  following  examples  will  be  considered  at  length  further  on, 

but  may  be  touched  on  briefly  by  way  of  illustration.  Let  us  consider  the  boundary 

conditions,  u  given  on  the  left  half  of  an  infinite  straight  line,  •^^-Sjgiven  on 

^  7  I 

the  right  half.  If  we  use  a  cartesian  coordinate  description,  we  may  say 


u(x,0)  =  f(T)        x<0^ 


-|^)^^^  =  g(x)    .^0^ 


on  y  axis 


This  is  a  two  part  problem.  However  in  polar  coordinates,  we  would  have 

u(r,Tr)  =  f(r) 

-^(r,0)  =  rg(r)       , 

and  here  the  probl«ii  has  been  broken  up  into  two  distinct  one  part  problems.  Thus 
the  same  state  of  affairs  is  a  mixed  boundary  value  problem  when  expressed  in  car- 
tesian coordinates,  while  in  polar  coordinates  it  is  a  simple  boundary  value  problem, 
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Consider  now  a  problem  referring  to  an  obstacle  sltunted  upon  the 
X  axie  from  x  =  0  to  x  =  1.  We  may  have  boundary  conditions  which,  when  ex- 
pressed in  cartesians,  will  read, 

u(x,0)  «  f(x)         -oo<x<0 
a^(x,0)  =  g(x)         0<  X  <1 

u(x,0)  «  h(x)  xy  1 

The  single   coordinate  line  y  =  0  is  here  divided  into  three  parts. 
However  in  polar  coordinates,   we  would  hnve 

ii(r,-n-)    -  f(T-)  all  r 

f|(r,0)   =  rg(-)  0<r<l 

u(r,0)   =  h(r)  r  >1 

The  first  of  these  equations  may  be  treated  independently.  The  remaining  pair 
constitute  a  two  p^^-t  mixed  boundary  value  problem. 

We  realize  also  that  the  amployment  of  elliptic  coordinates  would 
lead  to  a  simple  boundary  value  problem  here.   (To  thif  end  it  would  be  useful 
first  to  make  a  translation  of  origin  to  the  point  x  =  ^,  y  =  0,  and  a  change  of 
scale.)  Of  course,  in  order  to  employ  the  method  of  separation  of  variables 
under  these  clrcximst^rices  the  pwrtlal  diffe'-ential  equation  must  separate  in  the 
coordinate  system  considered,  i.e.,  there  mur>t  exist  product  solutions. 

U,  The  Method  of  Separation  of  Yariablas  Applied  to  Diffraction  by  a  Stral^i 
Edge, Regarded  as  a  Two  Part  Problam 

te   consider  aeain  the  boundary  value  problem  of  section  2.  Again  k  is 

taken  complex.  We  employ  the  same  arrangement  of  coordinate  axes,  and  ask  for  a 

solution  of  the  wave  eouatlon 
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In  the  form 


We  find 


V  =  F(x).G(y)  . 


1:114   *  5!M  ^  k2  .  0 
^M  B(y)  ^  ^   -  ° 

»"  2     r" 

and  hence  I—  must  be  a  constant,  say  a  ,  and  ii-  mst  be  another  constant 
F  G 

2 
say  -y3  ,  and,  also 

2     2    2 
-  a  -^4-  k   =  0 

or 

Solving  for  F  and  G,  we  find  that  solutions  possessing  the  character  of  a 
progressive  vmve  are 


1  «  X  :::  ly  k^  -  cx^y 


If  we  set  ex  =  k  cos  c^  ,  with  Cf^    real,  this  function  becomes 
glkr  cos  {o  -    Yqi    ^    (where  r  and  6  are  polar  coordinate  x  =  r  cos  6;  y  =  r  sin  9) 
l.e,,  a  plane  wsve  traveling  in  the  iD     direction.  We  are  not,  however,  re- 
stricted to  real  (0   ,  and  thus  of  may  have  any  range  we  pleare.  As  Ie  usual  in 
the  separation  of  variables  procedure  we  build  up  a  more  general  solution  in  the 
form 


■/ 


00 

f(c<  )  el"  '^l  11^  »^  -of  y  d  oc  . 


o()  ^ioc  xl  iV  k^  -o(^ 


-00 
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Thls  solution  represents  a  bundle  of  plane  iraves,  nnd  is  also 
essentially  the  starting  point  of  Sommerfeld^s  treatment  of  the  problem.  Let 
us  represent  the  solution  of  our  boundary  value  problem  as 

where  u  represents  the  scattered  wave.  We  shall  take  u  to  be  of  the  form 


if 


-oo 

where  the  upper  or  lower  sign  is  to  be  taken,  according  as  y  is  positive  or 
negative. 

Our  conditions  are 


(4.3)  ^^^  f? 


=  0  :  X  >  0 
7=0 


(b)     u  regular  in  the  exterior  of  the  screen. 

From  Eq.  (4.I)  and  (A. 2)  we  have 

^00 

U.A)  lim^      (u(x,€  )  -  u(x,-e  )^  =1    f(or)  e**^'^  d<X   . 

According  to  (4.. 3b)  u  is   continuous  in  the  exterior  of  the  screen,  and 
in  particular  on  the  negative  x  axis.  Hence  we  have 

00 
(4.5) 


I     f(cx  )  e^*  ^  dcx  =0  ;  x<  0 

-00 
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On  the  other  hend,  applying  (A. 3a)  to  (^^.1)  and  (4^.2)  we  find 


oo 
..   .  _    Ikx  cos  Oo 

Ik  sin  ©o  e 


1  r   f(o<')y  lc2  -0(2  e*  **  ^  dot   =  0;  X  >0 


-oo 


The  first  term  In  this  equation  Is  due  to  the  plane  wave  excltttlon. 
We  express  thla  term  also  as  a  superposition  of  terms  e      ,   and  Include  It 
under  the  Integral  sign.   (This  can  be  accomplished  by  the  Fourier  transform 
theorem.  We  only  need  the  representetlon  for  x  positive,  and  hence  may  let  this 
function  be  zero  for  x  negative.)  We  obtain,  using  the  fact  that  k  Is  complex 

(^•')   r  (1  !':  t'os  e„  *   I  r(«)/^^^^j  .'-^«  =0;  ..0  . 

^  "00  ^ 

Now  we  may  not  Infer  from  (A. 5)  and  (^.6)  that  the  respective  Integrands 
vanleh,  since  the  equations  hold  only  for  half  ranges.  We  seek  f(a)  fulfilling 
the  required  integral  conditions.  Now  in  order  to  evaluate  Integrals  one  often 
employs  the  method  of  contour  integration.  Let  us  then  regard  (4.. 5)  and  (4.6) 
as  contour  integrals,  completing  the  contour  in  each  case  by  an  Infinite  semi- 
circle in  of  plane,  where  o(  is  regarded  as  a  complex  variable.  These  semicircles 
will  give  no  contribution  provided  the  integrands  vanish  sufficiently  ra{;ldly. 
Now  along  such  a  semicircle  we  may  set  «  =  p  e    '       ,  where  ^Is  zero  on  the 
positive  real  axis.  Then 


e^*^^ 


=  e-'  P  ^^^  9 


Thus  the  exponential  factor  in  the  integrand  will  be  small  for  large  j>    , 
when  X  sintf?  Is  positive.  Hence,  for  x>0,  we  choose  a  semicircle  in  the  upper 
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half  plane,  sin  Op  >  0;  for  x<0,  we  chooBe  b   semicircle  In  the  lower  half  plane. 
We  now  oeek  en  f(c<  )  of  algebraic  growth  in  order  to  avoid  defeating  the  ex- 
ponential decay,  eo  that  our  semicircles  will  contribute  nothing  to  the  Integral. 
Now  we  may  view  the  Integral  in  (4-. 5)  as  a  contour  integral  where  the  contour 
encloses  the  lower  half  plane.  To  ensure  the  vanifhlng  of  the  Integral  we  seek  an 
f(  ex  )  regular  in  the  closed  lower  half  plane.   (Actually  weaker  conditions  on  the 
real  axle  will  suffice.)  Let  us  d«note  the  Integrand  In  (4-. 6)  ty  g(  V  ),  I.e.,  for 
ex  real,  we  define 


%     ^   sin  Oft         1    ,   X  I 2 2 
U.'')         g(«.  )  =  2 —   +  i  f(c/  )V  k^  -  cx^ 

o<  -  k  cos  e^^     '^ 

Elsewhere  g(oC  )  is  defined  by  analytical  continuation.  Then,  again,  we  can 
fulfill  (A..6)  by  requiring  g(o<  )  regular  on  the  upper  half  plane,  and  on  the 
boundaiy. 

Equation  (^.7)  is  manifestly  equivalent  to  the  transform  equation  In  the 
Integral  equation  method.  The  device  of  employing  the  monodromy  theorem,  analytic 
continuation  and  Liouville  theorem  appears  necessary  at  this  point  as  In  the 
Integral  equation  method.  The  overlapping  of  the  domain  of  regularity  is  here 
merely  on  the  real  axis,  but  this  is  sufficient  for  the  purpose. 

Coroparison  of  the  Present  Method  yith  that  of  Section  2:  Remarks  on  Tftilquenese 

.  One  may  note  in  this  heiirlstlc  procedure  the  manner  In  which  the  respective 
domains  of  regularity  ere  of  use  In  solving  the  problem.  Another  feature  of  In- 
terest may  be  noted.  The  conditions  at  infinity  are  rather  loosely  specified.  In 
parti cxilsr  a  point  often  glosned  over  in  the  discussion  of  the  problem  Is  that 
one  may  expect  a  reflected  wave  term  In  the  solution,  of  the  form  e^^''^'^°^'®o"''^^'^^o) 
in  a  suitable  region  of  the  half  plane  y<;0.  The  usual  uniqueness  requirement  is 


-3^ 

Ikr 
essentially  that  the  wave  u,  "be  of  the  form  ^-    ^(ft)^  I.e.,  that  It  decrease 

V  r 
with  increasing  r  In  this  manner.  Since  k  is  usually  taken  to  "l)e  real,  the 

reflected  plane  wave  term  has  not  this  hehavior.  But  with  k  =  k,  ♦  Ik^.  k 

i     2*   2 

positive,  one  may  require  instead  that  the  function  merely  vanish  at  infinity 
and  this  is  satisfied  by  the  reflected  wave. 

There  is,  however,  a  further  lack  of  uniqueness.  We  note  that  in  the 
present  heuristic  method  we  have  a  certain  latitude  in  discussing  the  asymptotic 
form  of  f(o^)  ar.d  g(<»<.).  When  we  reach  the  stage  in  which  the  entire  function  is 
represented,  in  the  various  half  planes,  by  the  left  and  ri^t  sides  of  a  certain 
equation,  we  can,  for  instance,  increase  the  exponent  of  OC  by  unity  in  the  asymp- 
totic forms  of  fit<  )   and  g(«sC  ),   (The  latter  was  designated  ^(«<)  in  Section  2). 
This  would  lead  us  to  an  entire  function  which  is  a  constant.  If  the  entire 

function  is  Z{pC  )  we  have 

«  ™/-^  A         2  i  k  sin  6 

(M)         f  (oc)  =  iM^  * 


i  Vk-oc  vlc+k  cos  0     (oc  -k  cos  6   )  vTo- ot 

0  0 


Here  the  second  term  on  the  right  is  the  previously  obtained  form  for  t_{oC)  ^ 
80  that  if  B(oC  )  is  not  zero  the  first  term  on  the  right  is  added.  We  note  that 
to  assume   the  growth  of  0(oC  )  in  this  manner  corresponds  to  a  non-integrable  be- 
havior of  ^(x)  near  the  origin.  In  the  Greens  function-integral  equation  procedure 
this  is  forbidden.  The  reason  is  that  there  the  asymptotic  form  of  the  transform 
is  deduced  from  the  behavior  of  the  function  at  the  origin.  This  in  turn 
is  limited  by  the  fact  that  the  intej^Tal  equation  oust  subsist,  and  hence  inte- 
grability  is  a  requirement.  The  Fourier  transform  faraula  also  required  Inte- 
grability,  but  this  may  be  overcome  by  deforming  the  contour  to  avoid  the 

singulnrity.  As  for  the  subsistence  of  the  integral  eountion,  it  may  be  ur.-^ed 
that  this  recrairement  is  logically  circular.  Certain  conditions  of  Integrability 
are  admittedly  required  in  order  to  v/rite  the  intet^al  equation,  but  on  the  other 
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hand  one  may  be  satisfied  to  retain  lack  of  Integra bility  while  rejecting  the 
obligation  to  wTite  the  integral  equation  without  et  least  some  modification.  It 
can  "be  shown  that  solutions  corresponding  to  non-integrability  of  f(:^  }  or  (8(7) 
at  the  origin  exist.  Let  us,  for  example,  consider  the  function 


V  = 


_1 ^_u 

~  ik  COB  9   ax 

where  u  is  our  previously  obtained  solution.  The  excitation  term  in  u  is 

gik(x  cos  60  + y  sin  %)  ^     ^^^   ^^^  definition  of  v  this  is  the  case  with  v  also. 

Also  -^-^  »  •: = — r-  -^: rr^  i  Hence  if  r—  approaches  zero  as  we  approach  the 

^y    ik  cos  9„   «3  X  Vyr         ^7     ^^ 


'o 


screen,  so  does  -4^  •  f^^   if  'tbe  scattered  part  of  u  vanishes  at  00 ,  so  does  the 
ay 

scattered  part  of  ▼.  Thus  v  appears  also  to  be  a  solution  to  the  problem,  and  if 
we  desire  uniqueness,  we  seem  forced  to  specify  the  singulerities  of  u  more  exactly. 
This  is  not  surprising,  since  the  uniqueness  theorem  states  that  a  wave  function 
for  an  exterior  problem  is  uniquely  specified  by  the  radietion  condition  and  by 
its  boundary  values  on  some  smooth  finite  boundary,  and  by  regularity  in  the  ex- 
terior of  the  obstacle.  In  cur  cases  the  boundaiy  is  not  finite,  the  obstacle 
has  a  sharp  point,  and  the  fxmction  is  not  regular  in  its  approach  to  that  point, 
With  regard  to  the  singularity  of  u  and  -^  near  the  point,  it  must  be  remarked     i 
that  they  affect  the  growth  of  f_(ec)  and  0^(oC),  and  hence  the  determination  of 
the  entire  function.  Furthermore  we  may  introduce  various  linear  combinations  of 
u  and  its  various  derivatives  with  respect  to  x,  corresponding  to  various  non 
integrable  behaviors  near  the  origin.  Thus,  to  insure  uniqueness  it  seems  suggested 
to  specify  the  behavior  at  the  origin  more  closely  in  advance. 

*  Conipare  also  Bouwkainp  l17j  i  where,  however,  the  relation  to  the  Viener-Hopf 
solution  is  not  In  rraestion. 
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In  illustration  of  the  above  considerations  consider  the  problen  with 
the  boundary  conditions 

u  =  0     at  7  =  0      x^O 

-f^  =  0   at  y  =  0    x>0 

u  continuous  In  the  exterior  of  the  x  axis, 
u  vanishes  like  an  outgoing  wave  at  infinity. 


A  solution  Is 


ikr 

.  cos       _   ,_, 

r        2      1        2 


u  =  S-J  .  cos  I  =  H^^\kr)  cos  | 


2ei^ 
and  this  has  the  Jxiinp  -^  •  •   on  the  positive  x  axis. 


V  X 

But,  other  solutions  are 

u  =  H   ,  (kr)  cos  (n-»- 1)  0  for  n  integral,  and  these  behave 

like  x~      for  y  «  0,  x  -*  -^  0.  In  actual  fact  these  can  be  derived  from  one 
another  by  differentiation.  We  have 

9%      ^  ^^       sin  e  9^ 

^^  =  cos  e  — ;:-  eo- 

and  using  the  well  known  recurrence  formulae  for  Hankel  functions 

^  ^  _  1  /~-     '-'         -i 

aT-^  K-i-%^1^ 

we  find  ^ 

Vi'Vi  --ST 
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It  appears  that  the  functions  u  may  be  added  In  various  comblnp.tlons  tc  our 
Solution  for  diffraction  by  a  straight  edge;  a  more  co.refvil  specification  of 
singularities  is  necessairy  to  exclude  them.  In  fact,  when  E(o«)  is  a  constant 
the  additional  term  in  (U,2)  corresponds  essentially  to  the  addition  of 


H|  (kr)  cos  I 


v;e  add  another  remark  about  the  significance  of  multiplying  f(«x)  by 
some  integral  power  of  ex  ,  and  the  relation  of  this  phenomenon  to  differentiation 
of  the  solution.  In  the  representations 


u  =  +  3-j    f(oc)e^  '-^^dof 

-co 


^'I 


/OO 
f(cx  )  e^'^''  d« 

-co 

we  may  proceed  operationally,  and  seek  to  differentiate  formally  with  respect  to 

»n 
X  under  the  integral  sign.  If  we  apply  (-^—j  ,  we  are  led  to  a  factor  (i  «)" 

multiplying  f((X  ).  Thus  the  occurence  of  terms  of  the  form   JP^  .  .  .  in  f^cx  ) 

Vk-  « 

may  be  regarded  as  symbol  1  zing  terms  likef^j— |  Hi(kr)  cos  -  in  the  solution  u. 
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5.  The  Method  of  Separation  of  Variables  Applied  to  Diffraction  by  an 
Infinite  Array  of  Semi-Infinite  Planes. 

In  this  section,  we  i»lsh  to  discuss  the  application  of  separation  of 
varitbles  In  Cartesian  coordinates  to  a  more  complicated  boundary.  The 
"obstacle"  in  question  Is  a  staggered  infinite  arrey  of  semi  infinite  parallel 
plates.  The  Green's  function  technique  introduced  by  Schwlnger  has  been  ap- 
plied to  this  problem  in  a  pair  of  papers  by  A.  E.  Reins  and  J.  F.  Carlson  ^^l^^. 
The  geometry  in  question  is  indicated  In  the  figure  beloti,  (figure  and  notation 
being  patterned  after  these  authors) . 


a^Tr<9<a<^ 

Direction    of 

Incident 

Wave 


T 

a 


a 


X 


\ 


^z 


a 


conducting    plate 


Fig.  4 


We  will  first  revlev  briefly  Carlson  rnd  Heine'  formulation  of  the 
problem  treated  in  their  second,  paper  [hB],  In  this  case,  the  Incident 
field  has  but  a  single  component  of  the  magnetic  vector  H,  and  this  component 
Is  the  component  perpendiculer  to  the  plane  of  the  paper.  This  is  then  true 
also  of  the  excited  field  due  to  the  presence  of  the  pletes.  The  components 
of  the  electric  field  E  lie  in  the  plane  of  the  paper,  >Jnd  if  the  x  com- 
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ponent  of  the  magnetic  field  is  denoted  by 


then  one  has. 


\(7»z)   -  )^(y,2) 


Ikl  =  - 


9   z 


»    ay 


idille    "yr  obeys  the  Helaholz  equation. 


yy 


r 


♦  k' 


zz 


y... 


0  . 


■^  is  subject  to  the  boundary  condition  ■^*-  »  0  on  the  plates,  since 

1   ^'^  ^ 

^z  "  7^  a    ^^  't^®  tangential  component  of  electric  field.  The  geometry  is 
lie  o  y 

certainly  not  a  likely  one  for  the  employment  of  separation  of  variables, 
Carlson  and  Heins  employ  the  Green's  function  technique,  using  the  free  space 
Green's  function  as  explained  in  section  2  of  the  present  report.  The  contour 
indicated  in  the  figure  below  is  introduced. 


contour 


Fig.  5 


dz' 
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In  an  array  of  m  plates,  and  then  the  radius  of  the  contour,  and  the  number 
of  plates  is  made  arbitrarily  large.  One  obtains  an  equation  for  j^  of  the 
type  encountered  in  section  2  of  the  present  report,  i.e., 

incident     m«oo  /  5y  L  J 

^  mq 

where  y'  =  a«  ,  q  -  d  cot«   and  U  (y,^)  -  e^^^^  ^^^  ®^  ^  °°^  ®^  .  The 

'inc 

quantity  Ij,(2*)  represents  the  current  on  the  plates.  Application  of  the 
boundary  conditions  then  leads  to  an  Infinite  set  of  simultaneous  integral 
equations  of  Fiener-Hopf  type,  i.e.,  each  possessing  a  difference  kernel.  Such 
a  set  would  not  readily  submit  to  the  technique  for  explicit  solution  de- 
•cribed  in  section  1  of  the  present  report.  The  fact  that  It  does  in  the  pre- 
sent case,  is  due  to  the  periodicity  of  the  structure.  Carlson  and  Heins  assune, 
i.e. 

I^  («  +  »q)  -  I^  («)  e^*™^**  '"'^   ^  *  "^  ^^"^  ®^   . 

With  this  assumption  a  solution  can  be  obtained,  and  verified  so  that  it  is  not 
merely  plausible,  but  correct.  Afer  a  change  of  the  variables  of  integration 
and  summition,  the  system  of  simultaneous  equations  is  reduced  to  a  single  in-  ' 
tegral  equation  with  a  rath--r  complicated  difference  kernel,  and  in  this  form 
the  now  customary  technique  can  be  applied,  as  Carlson  and  Heins  have  shown, 
despite  the  complexity  of  the  function  theoretic  analysis. 

I7e  shall  find,  as  is  only  proper,  that  a  similar  situation  obtains  with 
respect  to  separation  of  variables. 
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As  in  Section  4,  we  represent  the  perturb"^ tlon  due  to  each  plate  by 
means  of  a  sum  of  products  of  solutions,  choosing,  in  the  case  of  each  piste, 
an  origin  at  the  leading  edge  of  that  plate.  Thus  introducing  the  new 
variables 

z  =  z-mq  =  z-md  cot  <K  , 


ym  =  y-®^ 


we  write. 


where  the  signum  is  employed  in  order  to  permit  a  discontinuity  at  y  -  md  on  the 
plate.  Thus  we  have 


a^=  e^"*^^^  ^°^  ®o"^  y  sin  ^  ^  y2-   -^ 


m-   -co 


and 


^"^^  ik  Sin  9  e^^^"  "«  »o  ""  ^  '^^   ^^  7^  ^-^ 


sy 


On  each  plate  Te  must  h-jve     ^■'     =  0.      Hence  T/e  may   vrlte   for,    eay  the 


th     -    ^ 

n       piste. 


j1. 

5  y 


^0     =  ik  B 


.      -     ik(z  cos  e„*  nd  sin  G.  ) 

m  e„  e  0  o 


yn-  ° 


op  f    OO  y ^^ 

V  1       I  J  /»        4  ,/u2      /d2      r'a\      i-^v  z-ind  CO 

+       Z__     2TfT  ^^    '    *^^     ->^         V>^^    * 
m=  -OO            -/  -oo 


tc<  )+  i  I  n-m)  dVk^- j3^ 


for 


Zjjj  "7  0   ,    1 .  e .        z  ?  nq 


-^- 


We  have  one  such  equation  for  ever;  •»,  I.e.,  an  infinite  set  of  integral 
eqtifttionB  for  an  infinite  set  of  unknowns.  The  reason  we  may  succeed  here 
is  the  special  symmetry,  which  essentially  reduces  the  problem  to  one  in- 
volving merely  a  pair  of  the  plates,  i.e.,  an  ordinary  wave  guide  problem. 
Such  a  view  would  lend.  In  the  integral  equation  formulation,  to  dealing 
with  a  Green's  function  for  a  strip  shaped  region  and  this,  as  is  well  known. 
Is  obtained  by  an  image  method,  as  a  sum  of  an  infinite  number  of  Green's 
functions.  This  fact  is  reflected  in  the  separation  of  variables  procedure. 
1?e  recognize  that  there  will  be  a  discontinuity  of  the  function  on  each 
plate,  in  terms  of  which  the  current  density  may  be  represented.  Let  us 
therefore  form  the  jump  in  V  at  y  =  md.  1»e  have,  since  "^  continuous, 
and  'J^jj  continuous  for  yn  /  ^ 

r_,i         If    ^  f     \     i/3{z-md   cotoi  )  ^ 
J  -oo 

Hence,  (utilizing  the  fact  that  j^is  continuous  et  the  level  y  =  nd 
provided  x  <  md  cot cX  ,  i.e.,  off  the  plate),  we  have 

e-^""^  cotc<  f^(^)_  p      [y]  (^,^„d)  ^-i>9x'  ^, 

J   md  cot«t  ^ 

or,  If  we  use  s,  the  distance  from  the  leading  edge,  as  a  new  variable 
8  =  x'-nd  cot  cK  ,  we  have 

fn(.y3)~P   [y](s,n<i)  e^^^  ds  . 
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th 


Row  the  function  "Jr  her©  is  the  current  density  In  the  n   plate 
as  a  function  of  the  distance  from  the  leading  edge.  Ite  say  assuae  this  is  of 
the  form 


Ws.o) 


ikn(d  cote^   cos  0  +  d  sin  6) 


I.e.  of  the  same  magnitude  as  the  current  density  In  the  zero  plate  at  a  cor- 
responding distance  from  the  leading  edge,  emd  then  we  ba^e 

.     .  ,     ^      ikn  d  (cot«^cos  e.tsin  9) 

Consequently  we  have  for  ■^j-^  at  y  -  nd 

,  ^B,  :  aJ^  (z,nd)  =  0  -  ik  sin  0  e^*^^^  ^°«  V^  ^^-^  ®o) 
d  7 

00  

/  n=  -oo 

«'-oo  U 

with  q  =  d  eotof  as  before. 

Introducing  a  new  ynriable  z  =  z  -  nq.  I.e.  distance  measured  from  the  leading 

th 
edge  of  the  n   plate,  we  have, 

0  =  ik  sin  6  e^^^^  ''''^  V  ""^  ^^"  ^e^"^  *^"  ^  °°^  ®o 


y-OO 

-i-l   d>f  Oa)e^i^k^  ^^ikm(q  cosQ^^  sin  Q^)i\n-m\di^-^fi(n-m)oP 
/-oo    "  -oo 


for  z  >  0 


-u- 


or 


0=ik 


] 


with  Z'  -  q  COS  9,+  d  sin  9^ 

If  iM  write  B-n  =  p.  In  the  sum,  it  is  seen  not  to  depend  on  n.  The  sum 
may  be  written, 

p»  -oo 

and  representing  the  contribution  of  the  incident  field  by  a  Fourier  integral  for 
z  >  0  we  have. 


f'^  r  ~\    i>3  z 


Hence, defining 

we  seek,  as  explained  in  a  previous  section,  to  find  a  function  D(i3  )  which  is 
regular  in  an  apper  half  plane.  On  the  other  hand,  we  have  seen  that  f  (^ )  is 
proportional  to  the  transform  of  the  jump  in  ^'at  y  «  0.   Since  this  has  to  be 
zero  for  a  negative  we  h' ve  to  have  ^oifi)   regular  in  a  lower  hslf  ple.ne. 

The  stun,  K(y3)  is  readily  evaluated  by  tnenns  of  geometric  series  to 
yield 

cos  d-^  -y3^  -  COB  (ky>  -^q) 
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Hence,  our  eqtaations  ere  substantially  the  same  as  the  eouatlons  obtained  by 
Carlson  and  Heins  (equations  (3 •  A)  of  ref.  [l^.J)  except  for  a  misprint  which 
is  corrected  further  in  that  paper.   Our  fAj3)   correspond  to  one  half  their 
J(u; ) ;  our P  to  their  /^    , 

Thus  also  in  the  case  of  this  unlikely  looking  structure  one  is  led 
by  sepsrption  of  virifbles  to  the  eeme   position  as  by  aeans  of  the  Green *6 
function  plus  Green's  theorem,  as  it  was  desired  to  lllastrate.  Both  Bethods  con- 
verge on  this  point,  and  the  continuation  thereafter  is  the  same. 

In  the  H  plane  problem,  ^ere  the  boundary  condition  is  Y"  -  0  on  the 
conducting  plates  one  proceeds  in  a  similar  nay  to  the  same  end,  except  that 
the  perturbation  due  to  the  n   plate  is  now  written 


xoo 

'  n    2T1  i  I 

-'  -oo 


i  ct  (x-mq)  "f  1  Vk^  -a^  /  y-od  | 


f   d« 
ffi 


Here  the  omission  of  the  signum  (y-md)  makes  the  function  continuous,  while 
introducing  a  discontinuity  of  the  y-derivative.   This  must  be  removed  in  the 
exterior  of  the  plate,  but  is  allowed  to  remain  for  x-mq  ^  0. 

ThtLs,  by  either  method,  the  difficult  part  of  the  problem  lies  in  the 
necessary  function -theoretic  considerations,  and  has  been  solved  tor  Carlson 
and  Heins. 
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6,  the  iUthod  of  Separation  of  yarlatlea  Applied  to  Difftraetion  by  a  Strai<^t 
Idge  Begarded  a»  a  One  Part  Problem;  Solution  of  en.  Integral  Bqaatlon 

In  the  present  section  we  wish  to  consider  egain  the  botmdery  value 
problem  of  section  2.  In  sections  2  and  ^  this  problem  has  been  treated  by 
Wiener-Hopf  techniques.  We  heve  seen  that  the  basis  for  this  procedure  is 
the  fact  that  the  probleir.  is  being  regarded  as  in  cartesian  coordinetee,  whether 
explicitly  as  in  the  eigenfunction  treatment  of  section  4.,  or  implicitly  as  in 
the  Green's  function  integral  equation  treatment  of  section  2.  In  cartesian  co- 
ordinp.tes  the  problem  was  a  two  part  problem,  and  concomitantly  the  application 
of  the  Fourier  transform  to  the  Green's  f\inction  Integral  equation  required 
the  intervention  of  the  function  theoretic  TJiener-Hopf  techniques  for  its  suc- 
cess. (The  employment  of  the  Fourier  transform  is  of  course  regarded  as  the 
disguised  employment  of  carteeieji  coordinatee  and  constitutes  the  tool  for  de- 
composition of  the  Green's  function  into  the  exponentiel  cartesian  eigen- 
f\inctions) . 

On  the  other  hand  it  is  clear  that  in  the  polar  coordinates, 

x  =  rcose  o<e<2n 

y  «  r  sin  8  0  <  r  <  oo         , 

the  boundary  value  probleir  involved  is  a  one  p^rt  boundary  value  problem.     Ife 
restate  this  problem  accordingly 

(6.1)  ^    +   i      ^    +   i    ^    +   lc2u  =  0  :  Im.   k^O 

ar^         r       ar         ^.2     ^^^ 


(6.2)  u  =  Uq+-  u^ 


(6.3) 


_  lk(x  COB  9o  ■»•  y  Bin  Qq)  _  ikr  cos  (9  -  9^) 


(6.4) 


(6.5) 


a^i 

■*■ 

1 

r 

9^ 

^ 

ar2 

d 

r 

99 

r\»m 

= 

r\_ 

- 

2n 


0,  or 


U 


Ikr  Bin  Oft  e 


Ikr  cos  9, 


'0  =  0 


(6.6)      ui(r,9)  — >  0 


off  the  positive  x  axis 


oo 


(6.7)      u  continuous  except  on  the  positive  x  axis 


(6,8)      -~Q  continuous  evei^where  (except  possibly  the  origin), 


In  polar  coordinates  the  eigenf unctions  are 


Uv(r,e)  =  J^  (kr) 


\  cos  V  9  1 
(sin  ve  J 


where  V  may  be  complex.   Hence  a  more  general  solution  is 


1  oo 

•■■/ 

•^-1  oo 


J^  (kr) 


A(  V  )  COS  V  e  -J-  B(>»  )  sin  V  e 


dV   , 


and  one  may  expect  to  complete  the  solution  of  the  problem  ylth  this  repre- 
sentation in  a  clEcsical  manner,  not  Involving  function  theory.  To  do  this 
one  reoulres  the  following  transform  theorem  In  lieu  of  Fo\irier'£i  theorem: 


If 


(6.9) 


0(kr) 


^i  oo 
J  -i    oo 


V  ^  (V)   J^  (kr)  dV  , 


-^- 


then  _^ 

(6.10)  ^^^^=/,   -4^^  h/^)  (icr)  dr  . 

ThlB  theorem  has  been  devised  by  N.  H.  Lebedev,  and  employed  by  the 
latter  and  M.  J.  Kontorowltsch  [15]  in  the  solution  of  this  problem  In  the 
more  general  case  of  a  wedge  shaped  boundary,  (Actually  these  authors  con- 
sidered the  boundary  condition  u  =  0,  and,  furthermore,  applied  this  trans- 
formation directly  to  the  partial  differential  equation  after  certain  pre- 
paratory modifications  relating  to  IntegrabilityJ 

In  view  of  these  facts,  and  in  view  of  the  duality  between  the  Green's 
function  and  eigenfunction  procedures  we  may  consequently  expect  that  the  Green's 
function  integral  equation  (2.8)  of  the  problen  may  also  be  solved  by  the  em- 
ployment of  equations  (6.9)  and  (6.10)  without  recourse  to  f\jnctlon  theoretic 
methods.  The  Green's  function  Involved  is  a  constant  multiple  of  the  function* 


H^^^ky(x-xJ^-^  (y-yj2  ). 
o         o       o 

Introducing  the  quantities 


x^  =  r  cos  9- 
o    o     o 


yo""  ^o  ^^^  ®o 


this  function  may  be  represented  In  the  form 


oo 


H^^^k-/r2+  r^  -  2rroCos(e  -  9^,)  =  H   €„  H(^)(kr)  Jj^(kr^)co8  n(9  -  9^)  r  ?t^ 

n-o      ^ 

-°^     ,(2) 


=  IZ  ^n  ^^    H^^)J^{^t)   cos  n  (©  -  0  )   r<^  r 
n«o     n    ^     * 


with  6^  =  1,  ^„  =  2  for  n  >  0. 


»  fe   deal  with  the  function  H    henceforward  rather  than  H    here  in  order 

o  ^ 

to  accord  with  Ecs.  (6.9),  (6.10)  of  the  pr^;Eent  section. 
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To  form  the  kernel  of  the  Integrs.l  equation  we  need 


Ve   note  that 


a 

97 


y  ■  o 


-Id 

~  r   99 


0  =  0 


d    7  ijr. 


y  =  y^  =  0 


and  hence  our  kernel  Involves  the  function 


CO 


'^('^>-o>°^  -^  r:  ^„  n^  fl^^>(kx)  j^  (kx„) 

o  n=l 


(6.11) 


K(x,x„)oc  ^  IZ   «n  °^  H^^^l^o)  Jn(l«) 


o  n=l 


X  -?  x« 


x^  »o  ' 


where  XjX^  repl&ce  r,r  for  6  =  9o  =  0.  We  may  now  regard  this  sum  as  a  eeries 
of  residues  obtained  from  the  following  contour  integral. 

.1    00 

2 


(6.12) 


K(>^,Xo)  = 


const 

XX 


'/ 


X^  cot-iTA     H![2^(kxQ)  J^  (kx)  dA 


-i  oo 


The  contour  is  the  imaginary  axis.  It  may  be  deformed  into  the  contour  C 
illustrated  below. 


X-0 


Imaginary    axis 


X-plane 


• — *:  Real   axis 


X=i        X=2      x=3 


Fig.  6 
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to  obtain  the  pecond  form  of  (6.31)  when  x  <:  x  ,  The  integral  in  the  form 
(6.12)  is  sjnranetrlc  in  x  and  x^,  as  may  be  shown  by  repvesentlnp  the  Hankel 
function  In  tenne  of  Beseel  function  J>  and  J__^  .  Hence  when  x  >  Xq  ne  may 
Interchange  x  and  x^  In  (6,12)  and  obtain  the  first  form  of  (6.11)  upon  ex- 
panding In  residues. 

If  we  replace  the  kernel  in  the  Green *e  function  integrel  eauntlon  by 
this  integral  expression,  we  find 
(6.13) 


00  i  oo 

ikx  sin  So         "  =  const.  I   L—- ^U—S  /     i^cotTTA  Hi*^^  (kx„)J  (kx)d  A 


■^o      ""o  y.i   oo 


Interchanging  the  order  of  integration  and  defining 

oo 


o 
we  find 


f        [u  X-)j 
(6.U)  f(V)=/  ~^     H^(kXo)dXo 


ikx  cos  e^  /■  5  _ 

(6.15)       ik  sin  e©.  x*  e  =  const. I  ;\'^cotTr;^  f(A)J>  (kx)d,^ 


-i  oo 

anploylng  (6.9)  and  (6.10)  we  find 


r°°         ikx  cos  e^ 
(6.16)      AcottlA    f(  A)  =  const,  ik  sin  9^  I       ^^-^ " Hr^(kx)  d  A 


c 

I 


'Oj  X  "A 

o 


(2) 
To  evaluate  the  integral  on  the  right  hand  side  of   (6.]6)   we  represent  H  j.    (kx) 

in  terms  of  J,  (kx)  and  J_.  (kx).     We  ereploy  the  fonnul&6,[cf .    /?0],  p.  A05^  , 
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/(cob  bt  /  I  COB   (m  ?rcBln  —)  / 

J^   (at)  dt=    ^  ^  I    \ 

^     ^  ^8ln  btj  I  Bin  (;|  arcsln  |)  J 


These  fomrulae  are  valid  for  a  >  b:     Re  m>-1.     We  have^  -  ^  ^»   b  =  k  cos  9^, 
a  =  k,   so  that  a  >  b.     Also  we  may  deal  with  H\    (kx)  with  A  on  the  contovir  of 
integration  bo  that  Re  A  =  0,  whence  Reyl|  =  0.     We  find,   eaploylng   (6.16), 

-  —^    Bin  A  e 

(6.17)  Af(A)oC      e        2        -—^.^ 

whence,  employing  (6.9)  and  (6.10)  and  (6.14), 

A  <»  ITTA      g^         Q 

(6.18)  [u(xo)]oc/  e"       2       ^|_^    j  ^  (,,)  ,  ;^ 

^-1  oo 

This  expression  may  be  eveluated  by  expansion  In  residues.  The  result  Is  the 
series 

(6.19)  V'*^^)]^  TL     (-1)  J   .  (kx)  e       ^sin  (nvi)e^    . 

o         ""*•  ? 

This  may  be  expressed  in  the  customary  Fresnel  Integral. 

7,  An  aiectrostatlc  Problea  Relating  to  an  Infinite  Hibbon 
Introduction;  A  Two  ^e.rt  Problem  in  Polar  Coordinates. 

In  the  present  section  '>e  treat  the  m^.thematical  problem  arising  from 
the  perturbeticn  of  a  unifona  electrostatic  field  by  an  infinite  conducting 
ribbon.  We  take  the  z  axis  in  one  edge  of  the  ribbon,  and  the  y  axis  perpen- 
dicular to  the  ribbon.  The  ribbon  lies  in  the  xz  plane,  end  its  cross  section 
occupies  that  region  of  the  x  axis  which  lies  between  zero  and  one.  The  un- 
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pertarbed  ^leld  is  a  uniform  field-  and  the  lines  of  force  are  parallel  to 
the  X  axis.  The  following  figure  is  a  croBB  sectional  view  In  the  xy  plane, 
and  conditions  are  the  saiae  in  every  parallel  plane. 


y 

A 


unperturbed 
field 


0! 


ribbon 


■*-  X 


Fig.  7 

When  this  problem  is  considered  In  cartesian  coordinates,  it  turns 
out  to  be  8  3  pert  problem.  However  In  polar  coordinates  r,  0,  we  have 


X  =  r  cos  9 

y  »  r  sin  6 


0^9^  2TT 
0  <i  r  <  oo 


The  problem  •epare-tes  Into  an  Innocuous  one  part  problem,  plus  a  two  pert 

problem,  as  discussed  at  the  end  of  section  3.  Nor  the  el gen -functions  arising 

"♦•  V   ♦l  V  9 
in  polar  coordinates  for  Laplace's  equation  are  the  functions  r"    e 

Thus  a  solution  may  be  represented  as  a  sum  of  such  functions,  and  with  boundary 

conditions  at  fixed  values  of  9,  one  arrives  at  Mellln  transforms  rather  than 

the  related  Fourier  transforms.  The  latter,  on  the  other  hand,  occur  directly 

when  one  employs  cartesian  coordln^.tes,  for  the  product  solutions  are  e   \'^  "*  7) 

Thus  the  present  section  rill  Illustrate  the  transition  from  a  3  pert  to  a  two 

p?rt  problem,  and  the  solution  of  a  problem  by  the  Joint  employmect  of  function 

theoretic  techniques  and  sepprttion  of  varlsbles  in  a  case  vhere  transforms 

other  than  Fourier  trEnsfontis  are  involved. 
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If  the  problem  rere  set  up  as  an  integrel  eauation,  a  difference  kernel 
would  result,  but  the  range  of  integration  would  be  from  zero  to  one.   The 
Mellln  transform  would  be  suitable  for  effecting  a  solution  in  thle  case. 

Admittedly  this  problem  can  be  solved  by  conformal  mapping  in  a  simple 
■anner.  It  ie  given  here  as  an  llluptretlon  of  the  approach  of  the  present 
work. 

Solution 

Consider  a  uniform  electrostatic  field  In  the  horizontal  or  x  direction, 

and  the  effect  on  it  of  an  Infinite  plane  metsllic  ribbon  placed  rrith.   its  gen- 
erators at  right  angles  to  the  direction  of  the  field.  The  g-eometr>-  is  Indi- 
cated in  the  figure,  which  is  a  cross  sectional  view  in  the  x,  y,  plane.  The 
c  axis  is  taken  along  the  edge  of  the  ribbon. 

Fe  find  it  convenient  to  deal  with  the  stream  function  which  is  the 
conjugate  harmonic  function  to  the  potential  function  0.  If  F  ,  E  sre  the 
X  and  y  components  of  the  field  we  have 


(7.1) 


^   9x     9  y 


nr    ay      9  X 


and  of  course  0_,  +   (^™,     =  0     =  7^__     "♦■    /.„—   •      Since   the   tangential   componfri 
of  field  vanishes  on  the  conducting  ribbon,   we  have 


(7.2) 


,31 

3y 


=    0 


0  -c   X  <  1 


y»o 
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Also,  in  view  of  the  Bymmetry  about  the  x  axie,  EL  is  an  odd  function 
of  y.  Since  It  Is  continuous  on  the  free  part  of  the  x  axis,  it  must  venlsh 
there,  and  hence  so  must  -^*  .  Thus  ^  is  a  constant  on  each  free  part  of  the 
X  axis,  and  we  choose  this  constant  to  be  zero.  Thus  we  have 

1^=0   ,      fory  =  0     ,     x>l 
(7.3) 

y^-  0   ,      for  y-0     ,    -oo^x^O 

Alro,  the  field  at  infirity  should  be  une.ffected  by  the  ribbon,  and 

hence  E  =  1,  at  do,  if  we  consider  a  unit  incident  field,  or 
X    '      ' 

(7.^)  -^  (oo,  e)  =  1 


Fe  now  »Tlte   yj^  as  a  sxim  of  two  terms, 
(7.5)  /=7-|-H(x,y) 

where  H(x,y)   represents  the  effect  of  the  ribbon.     Then  we  have 


(7.6)  «xx  ^    V    =    ° 


(7.7)  H(x,0)     =0  ;         -oo-:.x<0 


(7.8)  H(x,0)     =0         t  1  <ix  <oo 


(7.9)  H  (x,0)     =     -1   ;  0<  X  cl 


(7.10)  II  (oo,  e)     =     0 


-55- 


We  nov>   separate  verlsbles  in  polar  coordinates  r,  9,  where  0  <:  0  •«^  2  TV 

X    =     r  COB  9 
y    ■     r  sin  0 


Equation  (7.6)  becomes 


and  solutions  are 


H   •*--Hf-LH   -0 


Hy   =  r   .• 


Superposing  a  combination  of  these,  we  write 

i  oo 


/I  oo 
r'''rA(v)co8V9i-B(Y)8lnV9J  dV 
-i  oo 

Then  by  (7.7)   we  have  for  all  r,  at  9  =Tr      ,    (i.e.,   the  negative  r  axle) 

i  oo 
l/2Tri       I  r""^  rA(v)   COS  VTTt  B(v)   sinvnl     dV       =0 

-1    00 

and  hence 

(7.12)  A(V)   cosVTT   \     B  (V)    sin  V  TT    =     0 

We  are  left  with  a  two  port  boundary  value  problem  on  the  positive 
X  axis,     re  hfive  by   (7.8) 

i  oo 


/I  oo 
-1    OO 


(7.13)  1/2X11    /  r    '^      A(V)   dV      =     0    ;         r  ^1 
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and  by  (7.9)  (since  on  the  positive  x  axis  ~-  =  -  -^  ) 


1  oo 
(7.U)  r-fl/2-nil  r"^     vB(V)dy     =0;     O^r-dl 


'/ 


-i  oo 

From  the  last  tuo  equations  we  may  not  infer  the  vsnlshing  of  the  integrends, 
of  course,  since  the  equations  are  not  valid  for  all  r.  Before  proceeding, 
let  us  modify  (7.14)  by  representing  the  quantity  r  in  the  range  0  <  r  <  1  , 
as  follows, 

1  oo 


mi  J 


-       -         ,  y.-^  d  V 

^  ■  2rti      I       "^    1  +  V 

-i   oo 


This  is  a  consequence  of  the  Mellln  transform  theorem  yrhich  states  that,  under 

certain  circvunstances,   if 

/or  +  i  oo 
r"^     f  (V)     dV 
CT   -i   oo 
then 


=/ 


oo 


f (V  )  =    /  f (r)  r^     ^     dr 


Tn  our  case  we  use  f(r)  =  r  (0  ^  r  <  1) ,  and  f(r)  =0,  r  >  1.  Having  thus  so 

to  speak  expressed  the  incident  field  in  eigenf unctions,  we  hsve 

r  1  00 

(7.U)'  2TrT  J  r""^  [VB(V  )+  l/ltV  1      dV=     0;      0<r-:l. 

-1   oo 
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We  now  employ  the  notion  of  contour  integration.  Te   may  close  the 
contour  of  integration  (here  the  imaginary  axis  is  the  V  plane)  by  means  of 
infinite  semi  circles,  provided  the  Integrand  vanishes  sufficiently  rapidly 
on  these.  We  writs 

V  =  Pe         f   -\^l  0  =  arg  V 

and  note  thnt  for  r  >  1,  \r      \      — ^  0  exponentially  with  p   — ^  oo,  provided 
cos  0  >  0,  while  for  r  ^  1,  we  require  cos  ^  <  0  to  produce  this  effect.  Thus 
we  may  attempt  to  close  the  contours  of  integration  In  (7.13)  and  (7.14)'  hy 
infinite  semi  circles  in  the  right  and  left  half  plane  respectively.  Then  if 
A("V)  is  at  Boat  of  algebraic  growth  In  the  right  half  plane,  and  V  B( V  ) 
behsves  likewise  in  the  right  half  plane,  the  contributions  thus  arising  will 
be  vani shingly  small,  and  ^-e  may  regard  the  integrals  in  (7.13)»  (7.14)'  as 
being  Integrals  around  closed  curves.  By   Cauchy's  theorem,  they  will  then 
vanish  as  required,  provided  the  integrands  are  regular  functions  in  the  re- 
spective half  planes.  Let  us  define  VB(V)i'  l/l^V  -  g^   (V)  ,  Then  we 
have:  A^  (  V)  regular  in  a  right  half*  plane, 

g  (v)  regular  in  a  left  half  plane, 
and  the  subscribed  4-  and  -  signs  are  intended  to  serve  i.s  reminders  of  these 
facts.  If  we  employ  equation  (7.12)  we  find 


le,^ 
•'  J 


(7.15)         g.  (V)  =  -VA^(Y)  ^^   -  ^    , 

if  there  exists  a  region  of  the  V  plane  in  which  all  these  quantities  are  regular. 
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The  procedure  is  now  similar  to  th^t  in  the  two  part  cartesian  case. 
Before  continuing  we  digress  to  establish  certain  results  which  will  be  neces- 
sary, We  recall  that  the  study  of  regions  of  regularity  and  the  study  of 
growth  at  infinity  in  the  V  plane,  plays  a  deeislTe  role  In  the  method;  we 
wish  to  indicate  the  form  these  investigations  adopt  here.  From  the  Mellin 
theorem  we  find,  (since  H(r,0)  =  0  for  1  <  r  <  oo) 

1 

(  v-1 

(7.16)  A^(V)  =  /    H(r,0)  r     dr   . 


If  we  assxuie  H(r,0)  '^  r~  *  near  the  origin,  then  the  integral  con- 
verges and  A^  (V  )  is  regular,  provided  Re  V  -  oc^>  0,  I.e., 

(7.17)  A^  (V)  regular  f or  Re  V  ^  W^  . 

On  the  other  hand  ^     (x,0)  =  "  3^  (r,0)  »  -  1  for  0  <  r  -i  1. 


Hence 


„V)    --/      :|fl  ,.,0,  .-1  d.  =/        II 


(7.13) 


or 


/I  00  /oo 

^''  ^^*f  If  ^"''^  =-(i/i4J   If- 
0        -^  1  1 


dr 


C-W  gJV)     --/        If    r^dr    . 
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Suppose  "g;r  -^  r     for  sufficiently  large  r.  Then,  as  far  as  convergence  near 
infinity  is  concerned,  we  have 

g  (v  )     regtilar  for  Re  V  ^  <•  -1 


We  have  supposed  that  ttie  Integrals  here  are  Integrable  at  the  other 
limit,  l.e,  at  r  =  1,  If  not,  ^re  cannot  of  course  regard  A.(V),  VB(v), 
g_(V)  as  Mellin  transforms  (unless  possibly  the  real  Integrals  In  (7.16)  and 
(7.19)  are  replaced  by  suitable  contour  Integrals),  but  a  solution  in  the  form 
(7.11)  may  exist  Independently  of  this  interpretation  ofA^("V),  'VB(v), 
g-(V). 

Iflth  regard  to  the  asjrmptotlc  form  of  A^  (V  ),  g  ( v)  we  make  the 
following  remarks.  The  integrals  of  (7.16),  (7.19)  contain  a  factor  r"^^  in 
the  integrand.  If  Re  V  =  0,  i.e.  if  Vis  pure  iaagin5.ry,  this  factor  os- 
cillates very  rapidly  (as  is  seen  from  the  form  r*    =  exp  V  (log  r)   )  for 
large  Y.  If  r<  1,  and  Re  V  >  0,  r*^   is  exponentially  damped,  and  this  is 
also  the  case  when  r>l  and  ReY  -<•  0.   In  any  case  at  oo  in  their  respective 
half  pl!ines  of  regularity  the  integrands  of  the  integrals  for  A.(V),  g  (y) 
approach  zero  rapidly  with  increasing  \V  I,  except  in  the  neighborhood  of 
r  =  1,  where  the  decay  or  oscillation  is  ineffective.   (It  may  be  rem3rked 
that  the  effect  of  the  rapid  oscillfition  in  the  first  mentioned  case  is  to  make 
alternate  infinitesimal  elements  of  the  integrand  cancel  each  other.)  Thus, 


for  the  growth  of  the  integrals  defining  A^.(V),  g_(V)»  the  neighborhood  of 
r  =  1  is  crucial.   If  we  suppose  that  fl(r,0),  ~—  (r,0)  behave  near  r  =  1 

-A  -4> 

like  (l-r)   ^  ,  (r-1)   "^  respectively,  we  can  say  that  asymptotically 


f       v-1         -A 

(7.20)         A   (V)~  j        r     (l-r)     dr  asV— »oo  in  the  right  half  plane, 
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oo 


>'rr-ir^2 


(7.21)      g  (V)'^  I    r  '^(r-1)  '  -^  dr       as  V -»  oo  in  the  left  half  plane. 


The  integrals  occurring  on  the  right  can  be  evaluated  explicitly  in  terras  of 
the  Beta  fiinction.  We  have, 

1 


B(p,q)  E  f      rP-1  (1-r)''-'   dr  =  ^jp)    Hq)   . 
J  0 


p(p+q) 


Thus 
(7.22) 


A(V) 


r(v)ra-A) 


If  we  make  the  change  of  variables  r  =  l/s  in  (21)  we  find 

1 


(7.23)    g 


._(V)^J   a  (1-,,)     d8=  FT^l ' 


To  evaluate  the  growt'i  of  these  expressions,  we  employ  Stirling's 
formula,  which  states  that 

(7.24)  r  {It  w)~  w"  e"'  V2V^ 

for  \wl  — >  00,  -Tr<  arg  w^TT  .  '^e  may  apply  this  formula  to  (7.22)  with 
w  =  V  -1  in  the  numerator,  w  =  V  -^3  ,  in  the  denominator,  and  obtain 


(7.25) 


-U. 


r      r(y) 


^.i^"-^  e-(^-^) 


^^'^   (1-1/V )''  >   e^ ^     v^  -  ' 
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Also  -V  Is  In  the  prop-r  range  of  arg  w,  when  V  is  In  the  left  half  plane,  and 
we  obtain,  from  (7.2A) 

(7.26)  p(.y)   --^  (-V)  2     • 

>^1-1 
Hence  a  growth  llk«  V      for  A  (V  ),  corresponds  to  the  behavior 

-A  A-i 

H(r,0) '^  (l-r)     near  r  =  1,  and  a  grovrth  like  (->> )      for  g_(V)  corresponds 
to  the  behavior  ^~   '^  (r-1)     near  r^ ,  These  inferences  as  to  correspondence 
apply  to  the  case  of  integrabllity  near  r  =  1,  i.e.,  1- >fl  >  0;  1- A  ^  0. 
Let  us  now  return  to  eouatlon  (7.15), 

g.^v;      VA^^V  )    gj^-^^  ^  14V 

In  this  equation,  A^(V)  is  regular  for  Re  V  >  <X  ,  g_(V  )  is  regular  for 

0^1  9  H  ~  °^ 

Re  V  <•  «=<  o  -  1»  where  H(r,0)  --v^  r"   ■'■  near  the  origin,  ^-   (x,0)  '^-'  r   ^  at 

Infinity.  The  half  plane  for  g_(V),  A^  (V),  will  overlap  when 

«1^  °^2-l   . 

In  our  case  we  specify  (X,  =  -  §,   ^  =2,  when  the  common  strip  of  regularity  is 
-  ^<Rey  <C  1,  enclosing  the  contour  of  integration  in  (7.11).   (Otherwise  one 
would  have  chosen  that  contour  to  lie  in  the  common  strip,  wherever  It  might  be.) 
Thus  we  refer  to  the  case  where  H(r,0) '>•' Vr"  near  the  origin,  and  as  a  result  this 
is  true  also  for  ^  .  If  we  want  the  same  behavior  at  the  other  end  of  the  ribbon 
we  have  to  choose  >^-  -  -  ^   also.  But  in  any  case  we  have  to  havey^^l,  J^2^^ 
for  integr^ble  solutions.  Solutions  corresponding  to  greater  values  of  y^  than 
chosen,  correspond  to  the  addition  of  various  derivatives  of  H  with  respect  to 
X.  These  hive  liigher  singularities  at  the  edges  of  the  ribbon,  but  do  not 
affect  the  boundary  conditions  on  the  plate  and  at  infinity  and  hence  a  solution  of 
the  problem  remains  so  when  these  are  added,  vt   least  as  far  as  equations  (7.5-7.10) 
are  concerned. 
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With  respect  to  the  choice  ol    =  2,  we  add  the  following  remark.  It  is 
well  known  that  in  the  case  of  a  uniform  exciting  field,  the  induced  field  be- 
haves at  Infinity  like  a  dipole  with  axis  oriented  in  the  direction  of  the  ex- 
citing field.  Thus,  at  infinity,  — -^  is  the  Induced  E  ,  and  hence  behaves 

9  y  x' 

like  l/x^  on  the  x  axis,  explaining  the  choice  of  cx  «  2, 

We  now  factor  the  coefficient  of  A^(  V)  into  factors  regular  in  the  right 
and  left  half  planes  respectively.  To  do  this  we  employ  the  Identities, 

vr(v)  =  r  (ifv ) 

cos  IT  y 

Substituting  in  (7.15)  we  have, 

Hulti plying  thru  by   p  jf  lyj  "e  find 
We  rewrite  the  last  term  as  follows; 


n 

r 


(1  -y)  *  i-vv    [  p(i  -V)    pTs)  Jl-^v  ^  p(?)(L  +  v)  » 


where  the  first  term  on  the  right  is  now  regular  for  Re  V  ^  ^  \*   since  the  pole  at 
V  =  -  1  has  been  nullified  by  the  zero  of  the  coefficient  at  the  point.  Then 
(7,23)  becomes 
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»   (V)    rjj^-v)  -  [  ^{^-y  I  _    ^  0/2)]  _i 


1__ 

V 


■»-^^^  p(f  tv)      r(2)      1  + 


The  left  hand  side  Is  regular  for  ReV'^  +  ^,  while  the  right  hand  side  is  regular 
for  Re  V  ?  -  J^.  Thus  each  side  defined  the  same  regular  function,  B(  V )  say,  in 
the  strip  -  ^  .^  Re  V<:  J,  and  the  continuation  of  E(  V  )  to  the  right  and  left  ia 
also  furnished  by  rtght  and  left  sides  of  (7,29).  Thus  E(V)  is  an  entire  function. 
Let  us  exaaine  its  growth. 

In  the  right  half  plane  we  have  A  ^(V  )'^  "^  ,  and  j^U  ^^  j  ^\^     V*  *  , 

in  accordance  with  (7.25).  Thus  the  growth  of  the  right  hand  side  is  like 

Examining  the  left  hand  side  we  find  g_(v  )  '^^  V   ^  '    '  j^  "V?  ^  V 
so  that  in  the  left  half  plane  the  representation  of  E(  V)  behnves  like  y^"^ 
Thus  the  entire  function  E(V)  is  asTmptotically  of  degree  less  than  one  if 
>3^4i  3/2,  J3^  <   5/2.   If  Integr^ibllity  of  H,  ~   is  assumed,  as  we  have  done, 
then  y^-t  /^2   *^®  both  less  than  1,  and  in  our  case  in  fact  y3    =  _  ^.  in  any 
of  these  cases  the  entire  function  is  then  certainly  a  constant  by  a  modification 
of  Liouville's  theorem,*  and  with  y3^  1,  this  constant  is  seen  to  be  zero.  Hence 

(7.30)         A^(v)=    r 0/2)   r;(»>v) ,  _j_.  To/a  siik±n 

(7.31,  ..(.,    =      ^-     f^-    I^-    -^t^ 


♦      Burlngton  and  Torrance:  Higher  Mathematics,   p.  602,  Theorem  17. 
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Hence  actuallj, 


A^(  V)  -^  V-  3/2 

g_  (V  )  ^  V-  i 


or 


A  =  -^ 


or  y9  =  +  i 


and  thus  H(r,0)  -^  V\^    near  r  =  1,  while  -^-^l  '^  .=?:-  near  r  =  1,  If  we  re- 


a_H| 


turn  to   {7.U)  we  find,  using  (7.30),    (7.32),    (7.12),   that 

i  oo 


H(...)  =  H(r.O,  =    IlUM  J  .-V    X|i^    ^i^'.JilL^    d  V 


(7.33) 


-i  oo 
i  oo 


=  r 


(^/2)       ( 

TTi         J 

-i  oo 


.-V 


VT 


1  tV)         P(i  -V)       C08VTT        ^V 


The  magnitude  of  the  ouantities  —^^^^17^^    '     '^'^'co^^-Tr  ^^       depends  on  that  of 

^i  V  p   -i  y  p 

h(V)  =   i^tr"   ivlT  »  ^^®''®  p  =TT  -  e.  Also,  in  our  problem  0^Q<2   TT . 
e    -  e~ 

This  quotient  is  clearly  bounded  for  fixed  Im  V,  as  ReV  — >  +  oo.  For  Increasing 

Im  V  =  T,  one  finds  the  following  asymptotic  forms  for  it: 


T  — }  ■>■  oo 

T  — ♦   -oo 

TT-  e  >  0 

h(v)-'  e"''^ 

h(v)'-  e*^ 

Tr-  e  ^  0 

h(v  )  '«  e     ^                 ' 

u^      X          T(2TT-  e) 
h(v)  '^e   ^ 

Thus  h(  V  )  is  bouiided  at  Infinity.  H(r,e)  can  thus  be  evaluated  by  expansion  In 
residues,  nnd  it  may  readily  be  verified  that  the  various  requirements  are  fulfilled. 
The  modification  for  a  ribbon  of  width  b,  and  a  field  of  strength  u^,  would  consist 
in  multiplying  by  u^^,  and  employing  r/b  Instead  of  r  in  the  integrals. 
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Summary  and  Remarks 

In  this  section  we  have  illustrated  the  transforaation  from  a  three  pprt 
problem  in  cartesians  to  a  two  part  problem  in  polars.  We  have  represented  a  solution 
of  a  pflrtlal  differential  enuatlon  as  a  sub  of  product  flmctions  by  separation  of 
variables,  and  h^ve  seen  that  the  situation  is  analogous  to  that  of  section  4.  The 
transform  theorem  in  question,  the  Mellin  Transformation,  is  inadequate,  but  suc- 
ceeds when  supplemented  by  function  theoretic  considerations  in  the  plane  of  the 
transform  variable. 

It  may  be  noted  that  the  present  problem  may  be  solved  by  alternative 
methods,  for  instance  conformal  mapping.  Also,  the  problem  may  be  reduced  still 
further  to  a  one  part  problem  by  the  employment  of  elliptic  cylinder  coordinates. 

Aside  from  its  value  as  an  illustration,  some  interest  may  be  attached  to 
the  present  treatment  when  it  is  regarded  as  a  model  for  ^  similar  treatment  of  the 
diffraction  of  a  plane  wave  by  an  infinite  ribbon.  In  this  case  the  equation  is 
^u'*lr-u  =  0.  A  solution  In  the  Mathleu  functions  arising  from  elliptic 
cylinder  coordinates  already  exists  in  the  literature.  This  corresponds  to  the 
treatment  as  a  one  part  problem.  The  use  of  polar  coordinates  should  lead  to  a 
solution  In  the  form  of  (7.11)  where  however  the  function  r~   Is  replaced  by  J  (kr) , 
To  complete  such  a  solution,  various  properties  of  J_^  (kr)  ere  seen  to  be  necessary, 
for  instance,  the  generalization  of  such  estimation  formulas  are  described  in  (7.21), 
(7.22)  and  the  preceding  paragraphs.  A  solution  has  not  been  obtained.  It  seems 
clear  that  a  modification  of  the  usual  function  theoretic  technique  is  necessary 
at  the  very  least. 
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8*  fhe  Pot«ntial  of  a  Ctharged  Dlao 

In  the  present  section  we  shall  Illustrate  our  thesis  by  considering  the 
problem  of  finding  the  elctrostatlc  potential  of  a  charged  disc. 

Let  the  disc  occupy  the  interior  of  the  circle  x^-^y^  =  1  in  the  plane 

2  =  0. 


Fig.  8 

The  problem  may  be  formxjlated  mathematically  as  follows.  We  week  a  f\m- 
ction  ▼(x,y,a)  such  that: 


(8.1) 


(8.2) 


▼   ••■  T   -^  ▼     SO 

XX   yy    88 


V  and  Vv  continuous  except  possibly  on  the  disc. 


(8.3) 


T  (oo)  S  0 


(8.A) 


T  s  ▼  ,  a  constant,  on  the  disc 
o'  * 


We  assume  axial  symmetry  about  the  z-axis  for  the  solution.  However,  we 
employ  spherical  polar  coordinates  r,  9,  J^,  where,  as  usual 


(8.5) 


X  =  r  sin  ©  cos  0 

y  «  r  sin  6  sin  fi 
z  =  r  cos  0 
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The  axlel  symmetry  Is  exprersed  by  considering  v  to  be  independent  of 
1^.     With  thlB  restriction  ft  solution  of  Laplace's  equation  in  the  fom 
^xii'*')  <L,(0)»   IJi   spherical  polar  coordinntes  is 

Y„(r,e)  =r-("*i)  P^Ccose) 

Superposing  solutions,  with  n-fl  =  V  ,  we  represent  ▼  in  the  fona 

<r-fri  oo 

(8.6)         '  '  2V1  /         '■"^  *^  ^>  V  -1^  JcoB  0))  d  V  . 

•^  0--  1  00 

Here  P    (cos  6)  is  the  Legendre  function.  When  the  index  Is  not  an  integer 
then  Py_i(cos  0)  will  possess  a  singularity  at  9  -TT  ,   and  P   .(-cos  9)  will 
possess  a  singulprlty  at  9  =  0,  The  use  of  P   .  ( |co8  8 1)  is  to  rcean  that  in  the 
upper  half  space  0  <  8  £  4^  we  use  P  ^  ^(cos  9)  while  in  the  lower  half  space 
5  ^  0  ^TT  »  we  use  P   ^(-coB  9).  While  P   .(Icos  9 J)  is  continuous  at  9  ="1  , 
since  the  argument  is  zero  in  both  cases,  this  is  not  the  case  with  its  deriva- 
tives, and  we  will  have  to  employ  special  care  to  avoid  dlscontlnultleE  of  V  ▼ 

off  the  disk  and  at  9  =  ^  ,  It  is  also  desirable  to  make  a  remnrk  about  the 

2 

quantity  <r  in  the  llirits  of  Integration  of  the  contov.-  Integral,  The  contour 
is  8   vertical  line  parallel  to  the  imaginary  axis  In  the  V  plane  at  e  station 
described  by  Ri  V  =  CT .  We  can  factor  out  r~^  from  the  Integral,  Since  ne  wish 
V  to  approach  zero  as  r  — >  00,  we  tentatively  set  P.e  O"  >  0,  to  ensure  this  be- 
havior subject  to  the  convergence  and  boundednens  in  6,  of  the  remaining  Integral. 

Clearly  ■  ^^  is  a  continuous  for  ©  =?  ,  and  thus  to  ensure  continuity 

d  v 

of  grad  V  we  hf.ve  merely  to  deal  with  -3^  .  Formally,  if  differentiation  under 

the  Integral  Is  permissible,  we  find 
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/(r-»-i  oo 
(T-  i  oo 


CT-t- 1  oo 


sLl 
do 


9^  J 


e 

2 

<J-  1  oo 


The  juBp  in  ^~   ,  denoted  by  j  "^    -tr  *  ^' 


[6  the  limit  of  the  difference  of 

2 
these  two  expressions  as  ^  — ^  0.  This  jump  is  to  be  zero  when  ire  are  not  on  the 


disc,  i.e.,  at  9  •  ^  and  r  >  1,  Hence  we  have; 


a-*"  i  00 

(8.7) 

Cjr-  i  oo 


a-*"  i  00 

~f  I         r-^A(V)  P'y_^  (0)  dV=0  :  r>l 


It  will  be  noted  that  the  r  in  this  expression  is  measured  in  the  plane 
z  =  o.   Hence  r*  =  x*-^.  y^  and  r  is  identical  with  the  radius  vector  of  cylindrical 
coordinates  r,  0,  z. 

We  have  now  formally  fulfilled  conditions  (8.1),  (8.2),  and  (8,3).  To  ftd- 
fill  (8.4.)  we  must  have,  by  (8,6) 


0-+  1 


C+  1  00 

r-V  AfM^  P 

V-1 


r-VA(v)P,.  ,(0)dV-Vo  =  0   :  r  ^  1. 


00 
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We  nov  digress  briefly  Into  the  theory  of  Legendre  fanctlone  In  order  to 
evaluate  P  -1(0),  P  ,(0),  We  shall  then  return  to  the  problem  of  determining 
A(v)  so  as  to  fulfill  conditions  (8.7)  and  (8.8). 

It  ic  known  that  P  y(cos  9)  may  be  represented  as  a  hypergeoraetrlc  function. 
Substituting V-1  forV  In  the  formula  for  this  representetlon  (cf.  [iS^  p.  311) 
we  have 

(8.9)  P^  _3^(coB  e)  ^  F(V  ,  1-V;  1:  1—^^^)       , 

The  recurelon  formula  for  Legendre  functions  is  (cf .  [l8j  p.  309) 

(8.10)  (z^-l)  P^  (z)  =    nz  P„  (z)  -  nP^  _^  (z) 

One  also  has  the  following  formula  due  to  Kummer   (cf .    (18 1  p.  298) 
(8  11)  F(2«     2/3-CX4     /9*i-i)     -      r  (*>*■•>  At   lA)    rL(l/gI 


Now  if  we  set  9  =  ^  in  (8.9)  we  find 


P^_3^(0)  =  F(V,1-V   ;1;  1/2)      . 

Setting  V=  2oC  ,1-V  =  2^   ,  we  find  <X  +^  +  l/2  =  1,    so  that  (8.11)   is  applicable. 
Using  r  (1/2)  ~  VtT  ,  we  find 

(8.12)  p,.  (0)  = ^-^ — 


V-l'      P(l/2  i-   V/2)  Td-  V/2) 
On  setting  z  =  cos  9*0,  corresponding  to  9  »  TT/2,   we  find  from  (8.30), 

P'v  (0)  -  VP^_^(0),  or  Py_i(0)  =(V-3)  Py.2(0) 
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or,  using  the  form  of  (8.12),  we  find 


(8.13)         P'    ,(0)  =z:       _        1      ^ 

V-1  P(V)p  (1^1.^) 


/tT"   (i/-i)     _   _  2-/Tr 


r(|)r(i-|) 


Here  the  formula   f  (1t^)  =    ^  r(5)  has  been  used,  P  denoting  ,   a«  usual,  the 
Gamma  function.     Hence  re  can  substitute  into  (8.7)  and  (8.8)  and  have  finally 
(8.7)'       ^O-^^oo 

2TT 


7)t       -Cr+i  oo 
■'  (T-  1  00 


.-V 


3^ 


''A(V)     p^^-.>d) 


dV  = 


0       :     r?  1         "V 


(8.8)'      ,^  +  ^  0° 


1 

2Tri 


^CT-  1 


-V 


0"-  1  oo 


A(v)  yrf 


r(^pr(i-p    ^ 


d  V  =  0:       r  <  1 
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Here  we  remark  that  we  have  replaced  the  constant    v      by  the  representetion 


'o         21^1 


^  0"+  i  oo 
^  <J-  1  oo 


r"^    V    dv 
o 


0<T^1 


with 


v/!'o'^-^-r°-^-^ 


civ     =     -^       (regular  ReV  >   0) 


To  fulfill   (8.7) *we  postulate 


(8.U) 


A(V)  n/TT 


r(V/2)P(l/2-V/2) 
and  to  fulfill  (8.8)'  we  postulate: 


Regular  Re.V  ^  CT 


(8.15) 


A(V)/X. 


r(|*^)r(i-|)     ^ 


-  -^  s  g(v) 


regulsr  Re  V  f  O" 
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together  with  proper  asymptotic  behavior  at  oo  In  the  N  plane  (i.e.,  algebrtiic 
growth  of  a  proper  order)  to  ensure  the  vanishing  of  the  integrals  of  these  Quan- 
tities on  sufficiently  large  circular  arcs  to  the  right  and  left  respectively  of 
tne  vertical  contour  of  Integration  In  the  plane.  We  note  that  a  contour  in  the 
right  half  plane  Is  suitable  f o  -  (8.7)'  since  r  is  greater  than  1  there.  The  re- 
verse is  true  for  (8.9) ♦.   It  nay  be  remarked  that  for  a  disc  of  radius  b,  one 

■Ight  replace  r  by  ^  in  (8.7)*  and  (8.8)*  so  that  the  coefficient  of  r~^  would 

b 

V 
be  multiplied  by  a  factor  b   ,  and  the  condition  r>  1  say  would  be  replaced  by 

r>b.  Rigorously  the  left  half  plane  considerations  night  requi^^e  special  treat- 

■ent,  but  of  course  one  nay  be  content  to  discover  an  answer  which  is  then  subject 

to  subsequent  verification. 

We  now  examine  (8,15).  Let  us  set 


(8.16) 

Then  we  have 
(8.17)  g 


My)ffr 


P(V)P(l_1t) 


(V) 


Jv)  =  -^      +   J  (V)  .    -TJ^l 


?^\^\) 


Multiplying  through  by 


ra-^j) 


T-    we  find,  at  Pe   V   =   0" 


P  (i  -^-) 


r  (1-  ^) 


(8.18) 


or 


(8.18) 


«.(^)ri^...!o 


I  2  2 


r  (1-  v>/?)  _  r  (1 


p(i/2-v/2)    TTi;^ 


fH 


)/     pa/2)   p  (1  i) 


2  2 
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In  (8.18)'  the  first  term  is  now  regular  for  Re  V  ^  cr ,    since  the  ap- 
parent pole  at  V  =  is  compensated  for  by  the  zero  of  the  bracket  factor.  Trans- 
posing we  find 

g  (V)   r,{Yi/2),:o\  riyi/^.) P ,(.ii1=  -Js.    n m ^ ^^wTj^ 

^-'"^^  p(l/2-v/2)*V    r(l/2-V/2)    P  (1/2)    V   p(l/2)     p7l72r?72) 

The  left  hand  side  of  this  equation  is  regular  for  Re  V  :^  CT  ,  and  the  right 

hand  side  for  Re  V  ^  (T  ,  provided  we  restrict  0"  still  further  b7  requiring  G"  4.    2. 

This  amounts  to  an  assumption  on  how  v  ->  0  when  r  — >  00  .  Then  an  entire  function 

E(V)  Is  defined  by  this  equation.  If  JzT  (  V ) ->^  T '^  with  y9<  |  ,  at  Infinity  in 

the  right  half  plane  and  ^(V  )  ^-^  Y'""°wlth  «  <  -  ^  ,  the  entire  function  la 

zero  by  Llcuvllle's  theorem,  using  suitable  Stirling's  theorem  estimates  for  the 

gamma  function  o.uotients.  Hence  we  find 

This  is  the  coefficient  of  r"^  in  eq.  (8.7)'.  By  use  of  (8.16)  A(  >»)  is 
determined.  We  now  wish  to  compare  this  method  and  result  with  one  due  tc  Titch- 
marsh  (Theory  of  the  Fourier  Integral,  Oxford  1937,  p.  334-  ff.).  Briefly  sketched, 
Titchmarsh's  procedure  is  as  follows.  Cylindrical  coordinates  y,  <3r,  z  are  intro- 
duced, and  then  v  satisfies  the  equation  (using  the  axial  ssmunetry) , 

with 

X  -   fcosCf  I     J  -    P  Bin  Cf ; 

The  Hankel  transformation  is  introduced,  one  defines, 


tu.,,  -/^ 


V(u,z)  =  1     f   w(f  ,z)   Jq  (j5  u)  d^>  ,       2  ;>  0  , 


_7 
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and  this  leads  to  the  equation 

V(u,z)  =  A(u)e   -►  B(u)  e 
and  B(u)  =  0,  for  Bmallness  at  infinity. 
Thus  actually  one  may  write 

00 


v(P,z)  =/  uA(u)  e""*  JqV^^  *'•  "  ^°* 

jThls  is,  in  other  words,  the  representation  obtainable  in  cylindrical 

coordinates.  If  z   Is  replaced  by  i a  J  one  has  a  representation  valid  for  all  z. J 

Friting  u  A(u)  =  f(u),  Titchmarsh  finds,  from  the  boundary  conditions, 

oo 
(     f(u)  Jo(j>u)  du  =  v^  :  0  <.f£\ 
o 

/oo 
u  f(u)  J^  (Pu)  du  =  0  ;  ^^1 
o 

In  this  position  of  ai'fairs  the  Mellin  transformation  theory  is  employed, 

and  using  the  convolutlonal  type  of  the  equations,  and  the  formula  for  the  Mellin 

transform  of  the  Bessel  function,  one  has,  where  F(s)  is  the  Mellin  transform  of 

f(u). 

k-t  i  00 

O^k  ^l 


1 

2irr 


rk  1 1  oo 

J  k  -1  oo  • 
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or,  with 


'<"  -=  -^^^}im-  X  (.) 


kM  oo 


1   I         I  (s/^)      -v      s-1 


i-k  +  i  oo 


At  this  stage  certain  ingenious  manipulations  are  undertaken.  The  first 
equation  Is  multiplied  by  P    and  integrated  from  P=Otop=lto  obtain 


fki-i   oo 

"^  i  k  -1  oo  1^^~^      "^  ds  =  ^-^.  .    Re.  u.  ^.  k 


Tr;r*v2r   '^^  °'  =  rr^;}  »     Re.  u.  ^  k  , 

similar  treatment  is  given  the  second  equation,  with  the  interval  1  4L  P  <    oo. 
Translations  of  the  contours  of  integration  and  the  employment  of  Cauchy's  Integral 
Theorem 


f  iU)) 


then  lead  to  the  result,  assuming  suitable  conditions  at  Infinity. 
One  finds. 


^o 


VT  (1-8) 

B-l 

Hence  the  coefficient  of  fi  in  the  second  integral  is 

^o       r  (1-  s/2) 

y]f  (1-8)    p  (l/2-s/2) 


-7^ 


If  xe  set  B  -  1  =  -V   ,  Re  s  =  k,  Re  V  =  or  we  find  0  <•  CTl.  1,  and 

the  coefficient  goes  into 


which  is  the  result  previously  obtained  (cf.  (8.7)',  (8.16)  (8.19)). 

From  our  previous  viewpoint  one  may  now  see  into  these  manipulations. 
The  introduction  of  the  Mellin  transformation  can  be  characterized  as  a  transition 
to  spherical  pol^r  coordinates  at  a  delayed  stage.  Since  if  r  is  the  spherical 
radius  vector  and  P  is  the  cylindrical  radius  vector, 

z  =  r  cos  0 
^  =  r  sin  0  , 

this  is  in  a  sense  a  tr.«>nsition  from  Cartesian  to  polar  coordin'ites  in  the  P  »  z 
plane,  for  the  eouation 

">/  +  p  7  +  ^zz  =  °    • 

We  see  that  the  |  functions  that  ST-ise  nre  the  reflection  of  the  Legendre 
functions  proper  to  the  spherical  polar  coordinates.  The  integration  with  l> 
and  the  use  of  Cauchy's  integral  formula  are  the  type  of  procedure  used  in  one 
proof  of  the  Mellin  inversion  formula,  and  may  be  viewed  as  an  alternative  to 
the  direct  employment  of  this  formula  involved  in  replaclni;  v^  ty 

fS\\.   oo 


.  _  -i  oo   ^ 
in  the  Interval  O^r^.1. 


-V 

^    d  V 


-%- 


Various  generalizations  are  poesible.   In  problems  where  Incident  fields 
are  involved  und  where  there  is  no  axial  symmetry,  one  might  represent  the  field 
in  the  form 


^    y'v      (r,e)  e^^f 
/  ..I  m 


in  spherical  polar  coordinates. 


Then  v      will  be  in  the  form 


TTj     r-^    Aji/)   Pj  (Icos  9|)   di/ 


''m  ■  2TT 


where  Pp  (cos  9)  is  an  associated  Legendre  function.  Application  of  the  boundary 
conditions  will  then  yield  eauatlons  like  (8.7)*  and  (8.3)',  but  with  more  compli- 
cated gamma  function  fftctors. 

The  following  is  another  extension.  Instead  of  r^.   disc  (or  aperture), 
one  could  consider  the  problem  presented  by  the  boundary  (0<r<l)  or  (r>l), 
and  ©  =  Qq,     This  is  a  hollow  truncated  cone,  (e.g.,  en  ice  cream  cone),  or  its 
complement.  For  Qq  =  ^  the  figure  degenerates  into  either  a  disc  or  aperture. 
In  this  case  the  formula  (8.6)  would  be  replaced  by 

cr+  1  oo 

r"'  A(y)  P     (cose)  dv       0^e^9- 

V  -1  *" 

(X-v  i  oo 

2^i  I         r"'  B(v)  P     (  -  cos  e)  d  V    .  e  ^.Q  <~   TT 
'   /  ff  -i  oo  '  "-'• 


/0-+ 1 
a-i 


vr-l--  I  .-V 
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The  treatment  Is  analogous  to  the  treatment  of  the  present  section,  but 
the  simple  formulae  for  P    (0)  are  no  longer  available.  Rather  one  has  to  deal 
with  P^  _i (c^e  9o).  The  crucial  features  are  then  (a)  the  expression  of  this 
function  as  an  infinite  product  in  the  i)   plane  fin  analogy  to  the  gamma  function 
arising  in  our  present  case  (9  »  -)L  (b)  the*  factorization  of  thla  product  into 
factors  regular  and  of  algebraic  growth  in  right  and  left  half  planes,  and  (c) 
the  determination  of  this  algebraic  growth.  This  work  forms  the  subject  of  a 
separate  report. 

If  a  siuilar  analysis  could  be  carried  through  with  Bessel  functions 
J^  (kr)  rather  than  p)ower6  of  r,  then  the  representation  of  a  solution  of  the 
scalar  w3ve  equation 

^at  1  oo 


A(v  )   ^'^  /L p^  (cos  e)  dv 

ar-i  00       V  5^     ^ 


would  lead  to  a  solution  for  the  acoustic  diffraction  by  the  above  boundaries, 
including  the  disc,  (cf  summary^  section  7  ). 


/ 
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